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Abstract 

Due  to  their  low  cost,  flexibility,  short  build-to-operations  cycle  and  capabilities, 
CubeSats  are  a  popular  approach  to  gain  access  to  space.  In  recent  years,  interest  in  larger 
scale  CubeSats,  such  as  3U  and  6U,  has  increased.  This  increased  interest  is  resulting  in 
new  challenges  including  generating  more  power  and  correspondingly,  radiating  more 
waste  heat  from  the  CubeSat. 

In  this  research,  a  modular  approach  was  developed  to  dissipate  heat  from  thermal 
energy  generating  components  in  a  CubeSat  in  order  to  prevent  excessive  temperatures 
and  meet  responsive  space  needs.  This  research  provides  the  details  of  the  technique, 
which  includes  physical  prototypes  and  implementation  methodology.  A  box  shape  heat 
sink  was  designed,  which  is  large  enough  to  cover  a  common  PC  104  circuit  board,  with 
the  purpose  of  dissipating  the  heat  on  the  circuit  board  through  thermal  straps  and  various 
size  radiators.  The  efficacy  of  the  design  was  demonstrated  using  both  computational  and 
experimental  methodologies. 

To  determine  the  heat  transfer  behavior  of  a  6U  CubeSat,  computational 
predictions  were  compared  against  experimental  results  collected  in  a  Thermal  Vacuum 
Chamber  (TVAC).  Subsequently,  the  validated  CubeSat  computational  model  was  used 
to  design  a  modular  heat  dissipater.  Validation  and  perfonnance  experiments  of  the  heat 
dissipater  were  conducted  using  prototypes  of  the  design.  All  computational  models  used 
in  this  research  were  validated  within  steady  state  with  Terror  <  4  °C.  Finally, 
computational  orbit  performance  predictions  of  the  heat  dissipater  were  made.  These 


IV 


predictions  revealed  that  the  temperature  of  a  heat  generating  custom  PC  104  circuit  board 
can  be  reduced  from  116.8  °C  to  12.2  °C  on  a  sample  hot  case  orbit  with  the  design 
created  in  this  research.  Alternatively,  it  is  possible  to  increase  the  power  of  the 
processors  on  the  circuit  board  more  than  400%  without  exceeding  the  temperature 
acquired  without  the  heat  dissipater  design.  Due  to  its  modular  approach,  this  design  can 
be  used  in  any  CubeSat. 
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MODULAR  HEAT  DISSIPATION  TECHNIQUE  FOR  A 

CUBESAT 

I.  Introduction 

Spacecraft  analyses  and  tests  are  extremely  important  for  exposing  potential 
failures  before  launch  and  ensuring  successful  on-orbit  operation.  Environmental  tests 
validate  that  space  systems  can  survive  in  the  harsh  space  environment  and  operate  even 
in  the  worst  cases  during  their  lifetime.  Thermal  testing,  a  part  of  environmental  testing, 
demonstrates  the  perfonnance  of  the  space  system  and  its  components  in  expected  and 
stressing  thermal  environments. 

The  Air  Force  Institute  of  Technology  (AFIT)  Center  for  Space  Research  and 
Assurance  (CSRA)  has  a  Thermal  Vacuum  Chamber  System  (TVAC)  with  a  solar 
simulator  sized  for  thennal  testing  small  satellites  such  as  CubeSats.  CubeSats  are  the 
new  generation  of  small  satellites,  developed  by  a  number  of  universities  and 
organizations  to  accelerate  small  and  low  cost  space  experiment  platforms  [Heidt,  2000]. 
Between  2003,  the  first  year  CubeSats  were  launched,  and  2012,  112  CubeSat-class 
missions  were  flown  [Swartwout,  2013].  In  addition  to  this,  197  CubeSat-class  spacecraft 
were  launched  in  2013  and  2014  [Swartwout,  2015].  These  data  show  that,  CubeSats  are 
gradually  providing  increased  access  to  space.  Today,  due  to  their  advantages,  CubeSats 
are  recognized  for  their  potential  utility  by  space  and  research  agencies  around  the  world 
including  United  States  (U.S.)  National  Aeronautics  and  Space  Administration  (NASA), 
U.S.  Department  of  Defense  (DoD),  National  Science  Foundation,  and  the  U.S.  National 
Reconnaissance  Office  [Asundi  and  Fitz-Coy,  2013]. 
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However,  Swartwout’s  analysis  on  CubeSats  shows  that  the  failure  rate  is  almost 
40%  of  all  attempted  launches.  In  addition  to  this,  launch  failures  are  also  20%  in  the 
years  between  2000  and  2012  and  the  failure  percentage  approaches  to  50%  in  university- 
led  missions  [Swartwout,  2013].  It  can  also  be  deduced  from  the  analysis  that  on-orbit 
failures  of  University-CubeSats  are  approximately  four  times  that  of  industry  mission 
failures.  Swartwout  concluded  that,  one  of  the  main  reasons  of  high  failure  rate  was  lack 
of  testing.  Through  analysis  and  testing  of  the  CubeSats  could  doubtless  reduce  these 
failure  rates. 

To  improve  thennal  analysis  and  validation  techniques  of  the  AFIT’s  future 
CubeSats,  a  computational  model  of  the  TVAC  environment  was  generated  by  Hatzung 
at  AFIT  in  2014  [Hatzung,  2014],  However,  this  model  was  validated  only  by  using  a 
10  inx  10  in.  aluminum  plate  because  of  time  constraints.  Since  the  TVAC  model  was 
used  in  this  research,  the  first  element  of  this  thesis  is  to  provide  a  greater  validation  for 
Hatzung’s  model  by  using  a  more  realistic  test  article,  a  6U  CubeSat  chassis,  for 
experiments  and  computational  simulations  with  a  steady  state  temperature  error  below 
4  °C  for  selected  set  points.  In  experiments,  dT/dt  <  0.5  °C  /  hr  was  determined  as 
thennal  equilibrium  and  steady  state  criterion. 

Another  fact  in  Swartwout’s  analysis  is  that  larger-scale  CubeSats  are  becoming 
more  popular,  because  of  the  need  for  larger  payloads.  Even  though  most  of  the  CubeSats 
have  been  lU-scale,  most  industry-built  CubeSats  have  been  3U-size  in  recent  years 
[Swartwout,  2013].  This  increased  interest  in  larger-scale  CubeSats  is  resulting  in  new 
challenges  including  generating  and  storing  more  power.  Correspondingly,  dissipating 
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the  heat  from  the  high  power  components  in  the  CubeSat  becomes  a  challenge  for 
thennal  control. 

Removing  the  waste  heat  from  the  components  is  critical,  because  it  provides  the 
opportunity  to  increase  on-board  processing,  improves  the  reliability  of  the  processor  and 
increases  the  expected  lifetime.  The  current  processors  on  satellites  are  operating  at  10% 
of  their  capability  because  of  insufficient  thermal  management  [Air  Force  Small 
Business,  2015].  Efficient  thennal  management  may  reduce  processor  temperatures  and 
increase  processing  capability  10  times  or  more. 

Research  Objective 

The  objective  of  this  research  is  to  develop  a  modular  approach  to  dissipate  heat 
from  thennal  energy  generating  sources  on  a  circuit  board  by  using  thennal  control 
hardware  in  a  CubeSat.  This  approach  can  reduce  time  spent  on  design,  building  and 
testing  of  the  CubeSat.  The  tasks  created  to  accomplish  this  objective  are  presented 
below. 

•  Task  1:  Provide  further  validation  for  the  previously  generated  TV  AC 
environment  computational  model  by  conducting  experiments  with  more  realistic  test 
articles,  because  it  was  validated  only  by  using  a  10  in  x  10  in.  aluminum  plate  as 
explained  previously.  Conduct  initial  experiments  with  a  6U  CubeSat  chassis.  All  models 
need  to  be  validated  within  steady  state  with  Terror  <  4  °C  compared  to  the  experimental 
measurements. 

•  Task  2:  Perform  computational  simulations  of  the  heat  transfer  behavior  on  a 
sample  6U  CubeSat  within  TVAC  and  compare  predictions  against  TVAC  experiment 
results  for  validation. 
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•  Task  3:  Create  and  analyze  a  modular  heat  dissipater  design  using  the  validated 
thennal  model  of  the  sample  CubeSat. 

•  Task  4:  Conduct  validation  and  perfonnance  experiments  of  the  heat  dissipater. 

•  Task  5:  Create  orbit  performance  predictions  of  the  heat  dissipater,  using  the 
validated  heat  dissipater  computational  model. 

This  heat  dissipation  technique  was  designed  to  prevent  excessive  temperatures 
for  the  heat-generating  component  with  a  plug-and-play  approach.  Due  to  its  modularity, 
this  design  can  be  used  in  any  1U  stack  in  a  CubeSat  and  reduce  the  design  to  launch 
period. 

In  the  remainder  of  this  thesis,  Chapter  2  provides  a  theoretical  basis  by 
presenting  background  to  heat  transfer,  thermal  environment  of  the  earth  orbit,  spacecraft 
thennal  design  and  spacecraft  thennal  analysis  and  testing.  Chapter  3  describes  the 
computational  and  experimental  methodologies  were  used.  Chapter  4  examines  the  data 
and  analyzes  the  results.  Finally,  Chapter  5  summarizes  the  results  and  offers 
recommendations  for  further  research. 
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II.  Literature  Review 


The  purpose  of  this  chapter  is  to  describe  the  theoretical  basis  for  the  research. 
First,  heat  transfer  and  the  modes  of  it:  conduction,  convection,  and  radiation,  are 
explained.  Second,  an  overview  of  thermal  environment  of  the  Earth  orbit  is  presented  to 
provide  a  basic  understanding  of  the  terms  direct  solar,  albedo,  and  Earth  infrared  (IR). 
Third,  spacecraft  thermal  design  process  is  expressed  and  some  commonly  used  thermal 
control  hardware  is  described.  Fourth,  spacecraft  thermal  analysis  is  introduced  and 
spacecraft  thermal  testing  is  briefly  expressed.  Finally,  an  overview  of  CubeSat  thermal 
design  and  analysis  research  in  the  last  decade  is  provided. 

Heat  Transfer 

As  it  is  well  known,  energy  can  be  neither  created  nor  destroyed.  It  can,  however, 
be  transferred  by  interactions  of  a  system  with  the  system’s  surroundings.  These 
interactions  are  called  heat  and  work.  While  thermodynamics,  the  science  of  energy, 
deals  with  the  end  states  of  the  process,  heat  transfer  provides  infonnation  concerning  the 
nature  of  the  interaction  or  time  rate  at  which  it  occurs.  A  general  definition  of  heat 
transfer  is  thennal  energy  in  transit  due  to  a  temperature  difference.  There  are  three 
modes  of  heat  transfer:  conduction,  convection,  and  radiation  [Ozisik,  1985]. 

Conduction 

Conduction  is  summarized  based  on  the  Incropera  et  al.  text  [Incropera  et  ah, 
2006].  This  heat  transfer  mode  is  transfer  of  energy  from  the  more  energetic  to  the  less 
energetic  particles  of  a  substance  due  to  interactions  between  particles.  It  occurs  in  a  solid 
or  a  stationary  fluid.  Rate  equations  quantify  heat  transfer  process  and  may  be  used  to 


5 


compute  the  amount  of  energy  being  transferred  per  unit  time.  Fourier’s  law  is  the  rate 
equation  for  heat  conduction  and  is  expressed  in  three  dimensions  as 

q  =  —kVT  =  —k(—i  +  —j  +  —k)  Equation  2.1 

dx  dy  d,z 

where: 

q  Heat  flux  (W/m2) 
k  Thermal  Conductivity  (W/mK) 

T  Temperature  (K) 

Vector  quantities  are  represented  as  bold  letters.  The  heat  flux  vector  q  is  the 

dT  dT  dT 

heat  transfer  rate  in  three  directions  per  unit  area.  — ,  — ,  —  are  the  temperature  gradients. 

dx  dy  dz 

The  parameter  k  is  a  thermophysical  material  property  referred  to  as  a  transport  property. 
For  a  prescribed  temperature  gradient,  increasing  thennal  conductivity  increases  the 
conduction  heat  flux.  The  minus  sign  indicates  that  heat  is  always  transferred  in  the 
direction  of  decreasing  temperature  gradient. 

In  case  of  a  thennal  energy  generation  within  the  system,  the  energy  source  term 
can  be  expressed  as 

Eg  —  q  dx  dy  dz  Equation  2.2 

where  q  is  the  rate  at  which  energy  is  generated  per  unit  volume  (W/m  ).  In  addition,  the 
energy  storage  term  is  represented  as 

Est  —  p  cp  dx  dy  dz  Equation  2.3 

where  p  cp  ^  dx  dy  dz  is  the  time  rate  change  of  the  sensible  (thermal)  energy  per  unit 

volume.  From  the  conservation  of  energy  law  and  the  first  law  of  thermodynamics, 
energy  is  always  conserved.  From  the  conservation  of  energy  requirement,  the  heat 
equation  can  be  obtained. 
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Equation  2.4 


a 

dx 


+  q  =  pcp 


dT 

dt 


Equation  2.4  is  the  general  fonn  of  the  heat  diffusion  equation  and  it  is  often  referred  to 
as  the  heat  equation.  This  equation  is  the  basic  tool  for  heat  conduction  analysis  in 

Cartesian  coordinates,  q  Is  the  energy  generation  term,  (p  cp  y-)  is  the  energy  storage 

tenn  as  explained  before.  The  spatial  and  temporal  temperature  distribution  T  can  be 
obtained  from  the  solution  of  Equation  2.4. 

Convection 

Convection  is  the  energy  transfer  between  a  surface  and  a  fluid  moving  over  the 
surface.  This  mechanism  of  heat  transfer  takes  place  because  of  the  motion  of  fluid 
relative  to  surface  [Ozisik,  1985].  The  equation 

q  —  h(Ts  —  Equation  2.5 

where: 

it  2 

q  Convective  heat  flux  (W/m  ) 
h  Convection  heat  transfer  coefficient  (W/m“  K) 

Ts  Surface  Temperature  (K) 

Too  Fluid  Temperature  (K) 

is  the  appropriate  equation  for  the  convection  and  referred  to  as  Newton ’s  law  of  cooling. 
Convection  heat  transfer  coefficient  depends  on  conditions  in  the  boundary  layer  such  as 
surface  geometry  and  thennophysical  properties  of  fluid  [Incropera  et  ah,  2006]. 
Convective  heat  flux  increases  with  increasing  difference  between  surface  and  fluid 
temperature. 
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Radiation 


Radiation  heat  transfer  mode  is  summarized  based  on  the  Incropera  et  al.  text  in 
this  chapter  [Incropera  et  ah,  2006].  Thermal  radiation  is  electromagnetic  radiation 
generated  by  thermal  motion  of  charged  particles  in  matter.  Radiation  is  emitted  by 
matter  at  a  nonzero  Kelvin  temperature.  Unlike  conduction  and  convection,  radiation 
energy  transfer  does  not  require  the  presence  of  a  material  medium.  In  fact,  vacuum  is  the 
most  efficient  medium  for  radiation. 

To  evaluate  the  radiation,  introducing  the  concept  of  the  blackbody  is  necessary. 
A  blackbody  is  an  idealized  physical  body  that  absorbs  all  incident  radiation,  regardless 
of  wavelength  and  direction.  Additionally,  no  surface  can  emit  more  energy  than  a 
blackbody  for  a  prescribed  wavelength  and  temperature.  The  Stefan-Boltzmann  law 
expresses  the  upper  limit  to  the  emissive  power  as 

Eb  —  o  Tg  Equation  2.6 

where  Eb  is  blackbody  emissive  power,  energy  is  released  per  unit  area  (W/m"),  Ts  is  the 
absolute  temperature  (K)  of  the  surface,  and  o  is  the  Stefan-Boltzmann  constant  (er  = 
5.67xl0~8  W/nr  ■  K4).  The  spectral  distribution  of  a  blackbody  is  represented  as 

2 hcl  Equation  2.7 

k,b(lT)  =  ^5[exp ((hc0)/AkT)  -  1] 


where  h  —  6.6  2  5  6x10  34  J-s  is  the  Planck  constant,  k  —  1.3805x10  23  J/K  is  the 
Boltzmann  constant,  c0  =  2.998x10s  m/s  is  the  speed  of  light  in  vacuum  and  T  is  the 
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absolute  temperature  of  the  blackbody  (K).  Blackbody  spectral  emissive  power  is  known 
as  Planck  distribution  and  represented  by  the  equation  below. 

E \,b  (A,  T)  —  nIA  b  (A,  T)  Equation  2.8 


Figure  1:  Spectral  Blackbody  Emissive  Power  (Evans  et  al.,  2014) 


A  plot  of  Equation  2.8  for  selected  temperatures  is  shown  in  Figure  1.  It  can  be 
noted  that  emitted  radiation  increases  with  increasing  temperature  at  any  wavelength  and 
the  sun  may  be  approximated  as  a  blackbody  at  5800  K  with  its  emission  peak  in  the 
visible  region.  It  can  also  be  seen  from  the  plot  that,  for  any  given  temperature,  Planck’s 


9 


function  has  a  peak  wavelength.  This  relationship  is  expressed  by  Wien 's  displacement 
law  as 

AmaxT  =  2897.8  pm-K  Equation  2.9 

and  it  can  be  deduced  that  with  increasing  temperature,  shorter  wavelengths  become 
more  prominent.  Whereas  a  blackbody  at  5800  K  has  its  peak  at  0.5  pm  in  the  visible 
region,  for  a  blackbody  at  300  K,  peak  emission  occurs  at  9.7  pm,  which  is  in  the  long 
wave  infrared  region. 

A  real  surface  emits  less  heat  flux  than  a  blackbody  at  the  same  temperature  and 
this  property  is  represented  by  the  equation 

E  —  £  oTg  Equation  2.10 

where  £  is  emissivity,  a  radiative  property  of  a  surface,  which  provides  a  measure  of  how 
efficiently  the  surface  emits  energy  relative  to  a  blackbody,  with  values  in  the  range 
0  <  £  <  1. 

Absorptivity  (a),  a  surface  radiative  property  of  a  material  is  used  to  evaluate  the 
rate  at  which  radiant  energy  is  absorbed  per  unit  surface  area.  This  rate  is  expressed  as 

Gabs  ~  a  G  Equation  2.11 

where  Gabs is  absorbed  incident  radiation  (W/m  ),  G  is  irradiation  (W/m“),  and  a  is 
surface  absorptivity  with  values  in  the  range  0  <  a  <  1.  Reflectivity  (p)  detennines  the 
fraction  of  the  incident  radiation  reflected  by  a  surface,  and  transmissivity  (r)  is  used  to 
determine  the  fraction  of  the  incident  radiation  transmitted  by  a  semi-transparent  surface. 
Each  surface  absorbs,  reflects  or  transmits  the  incident  energy.  This  can  be  represented 
mathematically  as 

aA  +  Pa  +  ta  —  1  Equation  2. 12 
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It  is  obvious  that,  if  the  surface  is  opaque,  it  either  absorbs  or  reflects  the  irradiation,  so 


ax  +  Px  —  1  Equation  2.13 

Additionally  Kirchhoff’s  law  states  that,  for  real  surfaces  at  a  particular  temperature, 
emission  is  equal  to  absorption  at  any  particular  wavelength,  which  can  be  expressed  as 

=  ax  Equation  2.14 


The  thermal  laws,  stated  in  this  chapter,  are  used  by  engineers  to  control  the 
temperatures  in  spacecraft.  Direct  solar  energy,  which  is  mostly  in  the  visible  spectrum,  is 
the  most  significant  heat  source  on  most  spacecraft  in  Earth  Orbit  [Gilmore,  2002]. 
Spacecraft’s  outer  surfaces  generally  emit  in  infrared  region  due  to  their  temperatures  as 
expressed  by  Wien’s  displacement  law.  Thus,  absorptance  in  the  visible  spectrum  and 
emittance  in  the  IR  region  determine  the  temperature  of  the  spacecraft.  Designers  select 
materials,  surface  finishes  and  coatings  of  exterior  surfaces  of  the  spacecraft  with  respect 
to  solar  absorptivity  and  IR  emissivity  values  to  set  the  desired  temperature  of  the 
spacecraft. 

Thermal  absorption  and  emission  of  all  inner  and  outer  surfaces  of  the  spacecraft 
can  be  selected  for  thermal  management.  IR  emittance  and  solar  absorptance  are  the  two 
main  surface  properties  for  surface  finishes  as  explained  above.  Since  the  main  heat 
sources  in  the  space  environment,  discussed  below,  are  mainly  in  the  visible  region  of  the 
spectrum,  materials  on  the  outer  faces  with  high  solar  absorptivity  increase  the 
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temperature  of  the  spacecraft.  On  the  other  hand,  outer  surfaces  with  high  IR  emissivity 
radiate  more  heat  to  deep  space  and  decrease  the  temperature  of  the  spacecraft  and  vice 
versa.  Thus,  the  designer  can  control  the  temperature  using  the  advantage  of  the  optical 
property  differences  of  surface  finishes  between  visible  and  IR  regions. 

Thermal  Environment  of  the  Earth  Orbit 

Environmental  heating  plays  a  major  role  in  spacecraft  thennal  control  and  the 
principal  forms  of  environmental  heating  in  Earth  orbit  are  direct  sunlight,  sunlight 
reflected  from  the  Earth  (albedo),  and  IR  energy  emitted  from  the  Earth  [Gilmore,  2002]. 
For  orbits  with  an  altitude  below  180  km,  free  molecular  heating  is  also  a  significant  form 
of  environmental  heating  which  is  a  result  of  bombardment  of  the  vehicle  by  individual 
molecules  in  the  outer  reaches  of  atmosphere.  Because  of  the  fact  that  free  molecular 
heating  affects  the  space  systems  especially  during  launch  ascent,  most  of  the  CubeSats 
are  assumed  to  be  thermally  safe  in  the  launch  vehicles  due  to  the  launch  vehicle  shroud. 
An  overview  of  direct  solar,  albedo,  and  Earth  IR  environmental  heating  forms  follows. 

Direct  Solar 

Direct  solar  is  usually  the  most  significant  heat  source  in  spacecraft  thermal 
control.  The  sun  is  a  very  stable  energy  source,  which  emits  the  radiation  remaining 
constant  within  a  fraction  of  1%  at  all  times  in  spite  of  the  1 1-year  solar  cycle.  However, 
Earth  has  an  elliptical  orbit.  The  intensity  of  sunlight  reaching  Earth  varies  between  1322 
and  1414  W/m  (approximately  +3.5%),  depending  on  Earth’s  distance  from  the  sun. 
Earth’s  average  distance  from  the  sun  is  1  AU.  The  solar  constant  is  the  intensity  of 
sunlight  at  1  AU  and  equals  1367  (W/m2)  [Gilmore,  2002]. 
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Albedo 


Albedo  is  the  fraction  of  the  incident  sunlight  reflected  back  to  space  by  Earth. 
Usually,  reflectivity  is  greater  over  continental  regions  than  oceanic  regions  and  varies 
with  local  solar  elevation  angles,  cloud  coverage,  or  ice  coverage.  For  a  spacecraft  in 
LEO,  the  orbital  average  albedo  varies  from  about  24  to  42%  of  direct  solar,  depending 
on  orbital  inclination  [Birur  et  al.,  2014]. 

Earth  IR 

Earth  reflects  some  sunlight  as  albedo  and  rest  of  the  sunlight  is  absorbed  by  the 
Earth  and  eventually  reemitted  as  IR  energy.  IR  energy  emitted  by  Earth  varies  with 
factors  such  as  local  temperature  of  Earth’s  surface  and  the  amount  of  cloud  cover.  For  a 
spacecraft  in  LEO,  the  orbital  average  Earth  IR  varies  somewhat  depending  on  orbital 
inclination,  but  is  generally  about  251  W/m2  [Birur  et  al.,  2014].  It  must  be  noted  that, 
since  albedo  is  the  reflection  of  sunlight,  it  is  mostly  in  the  visible  spectrum.  However, 
the  IR  energy  emitted  by  Earth  is  approximately  at  the  same  wavelength  with  the  energy 
emitted  by  spacecraft  [Gilmore,  2002],  For  this  reason,  Earth  IR  loads  incident  on  a 
spacecraft  cannot  be  reflected  away  radiators  with  special  coatings,  because  the  same 
coating  would  prevent  to  radiate  waste  heat  away  from  the  spacecraft. 

Spacecraft  Thermal  Design 

Spacecraft  thermal  design  is  a  process  to  control  the  thennal  energy  flowing 
throughout  the  spacecraft  system  in  such  a  way  that  thermal  requirements  are  met  during 
the  entire  life  of  the  mission.  Furthermore,  thermal  designers  should  make  the  optimum 
use  of  resources  such  as  mass,  power,  cost,  and  labor  [Birur  et  al.,  2014].  A  thennal 


design  process  starts  in  the  early  conceptual  phase  of  the  design  and  proceeds  through 
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various  stages  such  as  simple  heat  balance  calculations,  preliminary  designs,  more 
detailed  computer  analysis,  detailed  designs,  integration  of  the  hardware  and  finally 
thennal  testing  [Birur  et  ah,  2014].  The  main  purpose  of  all  efforts  is  to  set  and  maintain 
all  of  the  components  within  the  allowable  temperature  range  during  the  mission  lifetime. 
Establishing  a  thennal  design  for  a  spacecraft  is  usually  a  two-part  process.  The  first  step 
is  select  a  thermal  design  for  the  body,  or  basic  enclosures,  which  serves  as  a  thermal  sink 
for  all  internal  components,  and  the  second  step  is  to  select  thennal  designs  for  various 
components  located  both  within  and  outside  the  body  [Gilmore,  2002]. 

Besides  spacecraft  thermal  design  challenges,  CubeSat  thermal  design  presents 
additional  challenges  such  as  less  power,  less  volume  for  thermal  control  hardware,  and 
smaller  surface  area  compared  to  larger  spacecraft.  Although  a  wide  range  of  thermal 
control  hardware  and  techniques  is  available,  to  minimize  weight,  cost,  and  test 
complexity  and  due  to  power  and  volume  restrictions,  usually  active  thermal  control 
components  are  avoided  and  the  thermal  control  system  for  CubeSats  is  kept  as  simple  as 
possible.  To  provide  a  basic  understanding  of  spacecraft  thermal  design,  this  chapter 
introduces  the  common  thermal  control  hardware. 

Thermal  Control  Hardware 

A  wide  range  of  hardware  can  be  used  to  design  the  thermal  subsystem  and  meet 
thennal  requirements.  The  most  common  thermal  hardware  components  used  on  any 
spacecraft  are  multilayer  insulation  (MLI)  blankets,  radiators,  coatings,  and  heaters  [Birur 
et  al„  2014], 

MLI  blankets  are  used  on  the  outer  surfaces  to  prevent  both  excessive  heat  loss 
from  a  component  and  excessive  heating  from  environmental  fluxes,  rocket  plumes,  and 
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other  sources  [Gilmore,  2002],  MLI  blankets  are  typically  made  of  several  0.25  mil  thick 
metalized  sheet  layers  separated  by  a  mesh  made  of  a  low  thermal  conducting  material 
[Birur  et  ah,  2014]. 

Radiators  reject  waste  heat  to  space  by  IR  radiation,  and  occur  in  several  different 
forms  such  as  spacecraft  structural  panels,  flat-plate  radiators  mounted  to  the  side  of  the 
spacecraft,  and  panels  deployed  after  the  spacecraft  is  on  orbit  [Gilmore,  2002]. 

Heaters  are  active  thermal  control  hardware,  which  protect  components  from  cold 
temperatures  and  keep  them  within  allowable  temperature  limits.  Heaters  are  required 
when  passive  thermal  techniques  are  not  sufficient.  The  most  common  type  of  heater 
used  on  spacecraft  is  the  patch  heater,  which  consists  of  an  electrical  resistance  element 
sandwiched  between  two  sheets  of  flexible  electrically  insulating  material  such  as  Kapton 
[Wertz  et  ah,  2011],  The  patch  heaters  used  on  the  spacecraft  have  some  sort  of  control  or 
switch.  The  most  common  control  type  is  the  thermostatic  on-off  control  using  a 
bimetallic  mechanical  thermostat,  which  opens  or  closes  the  heater  circuit  at  a  preset 
temperature  [Birur  et  ah,  2014]. 

Surface  finishes  are  one  of  the  most  common  thennal  control  components. 
Various  materials  with  specific  wavelength-dependent  absorptivity  and  emissivity 
properties  are  used  as  thermal  control  coatings  with  different  purposes.  IR  emittance  and 
solar  absorptance  are  the  two  main  surface  properties  for  surface  finishes  as  discussed  on 
page  1 1 .  In  the  Gilmore  text  considering  these  properties,  thermal  control  surfaces  are 
fallen  into  four  basic  categories  with  respect  to  their  surface  properties:  solar  reflector, 
solar  absorber,  flat  reflector,  and  flat  absorber  [Gilmore,  2002].  Solar  reflectors  have  a 
very  low  a/s  ratio  (solar  absorptivity  /  IR  emissivity)  such  as  white  paint  and  second- 
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surface  mirrors.  On  the  contrary,  solar  absorbers  have  a  high  a/ £  ratio  such  as  polished 
aluminum  alloys  and  polished  gold.  Flat  reflectors  reflect  energy  throughout  the  spectral 
range  in  both  the  solar  and  IR  regions.  Aluminum  paint  and  polished  silver  are  examples 
for  this  type  of  surface.  Flat  absorbers  absorb  throughout  the  spectral  range  such  as  black 
paint. 

Spacecraft  Thermal  Analysis 

The  thermal  design  process  is  a  combination  of  design  selection  and  supporting 
analysis  as  emphasized  before.  Detailed  thermal  analysis  is  required  to  verify  and  refine 
the  selected  thermal  design.  The  thennal  engineer  starts  the  thennal  management  process 
with  establishing  thermal  design  requirements,  obtaining  spacecraft  design  and  details, 
determining  heat  dissipation,  and  identifying  environmental  and  boundary  conditions 
[Gilmore,  2002].  The  first  analysis  can  consist  of  simple  calculations  such  as  computing 
orbit  average  temperatures  and  power  dissipation  in  the  spacecraft.  However,  as  the 
design  progresses  and  becomes  more  detailed,  computer  simulations  are  needed  for 
analysis.  The  thennal  engineer  needs  to  construct  a  detailed  thermal  model  of  the 
spacecraft  to  predict  all  subsystem  temperatures  under  all  expected  enviromnental 
conditions.  A  description  of  a  typical  thennal  model  is  included  in  this  chapter  to  provide 
a  basic  understanding  of  thennal  software  calculations. 

Thermal  Model 

The  thermal  model  development  is  introduced  based  on  the  text  Spacecraft 

Thermal  Control  Handbook  by  Gilmore  [Gilmore,  2002],  A  detailed  thennal  analysis 

involves  constructing  a  geometric  mathematical  model  (GMM)  and  a  thennal 

mathematical  model  (TMM),  which  serve  different  purposes.  GMM  is  a  mathematical 
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representation  of  the  physical  surfaces  and  is  used  to  calculate  radiation  couplings 
between  all  surfaces  in  the  model,  as  well  as  heating  rates  to  each  surface  from  external 
flux  sources  such  as  direct  solar,  albedo,  and  Earth  IR.  In  this  research,  Autodesk’s 
Thermal  Desktop  (TD)  is  used  to  construct  the  GMM  of  the  6U  CubeSat.  TMM  is 
directly  analogous  to  an  electrical  resistance-capacitance  network  and  is  used  to  predict 
instrument  temperatures.  The  most  common  TMM  code  is  Systems  Improved  Numerical 
Differencing  Analyzer  (SIND A),  and  as  discussed  in  the  next  section,  it  is  also  used  as  a 
thennal  analyzer  by  thermal  desktop. 

Thermal  analysis  software  requires  the  user  to  define  a  thermal  network  system 
analogous  to  an  electrical  circuit.  To  achieve  this  requirement,  the  user  subdivides  the 
thennal  system  into  finite  subvolumes  called  nodes.  The  central  nodal  point  of  each 
subvolume  is  considered  to  represent  the  concentrated  thennal  properties  of  that  node. 
Each  node  has  two  thermal-network  elements:  a  temperature  that  is  analogous  to 
electrical  potential,  and  a  thermal  mass  that  is  analogous  to  electrical  capacitance. 
Thermal  analyzers,  such  as  SINDA,  usually  use  three  types  of  nodes  to  define  a  thermal 
network: 

•  Diffusion  node  (finite  thermal  mass)  has  a  finite  capacitance,  therefore, 
store  and  release  energy  over  time.  Thus,  the  temperature  of  the  material 
can  change  because  of  heat  flow  into  or  out  of  the  nodes. 

•  Arithmetic  node  (zero  thermal  mass)  is  a  physically  unreal  quantity  but  its 
effective  use  with  numerical  solutions  can  often  be  helpful  in  interpreting 
results  in  such  applications  as  surface  temperatures,  bondline 
temperatures,  and  node-coupling  temperatures 
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•  Boundary  node  (infinite  thermal  mass)  is  used  to  represent  a  boundary 
whose  temperature  is  set  and  will  not  change  such  as  deep  space. 

The  placement  of  the  nodes  depends  on  some  factors  such  as  the  points  where 
temperatures  are  desired,  the  expected  temperature  distribution,  and  ease  of  computation. 
Usually,  a  minimum  number  of  nodes  is  desired  to  verify  design.  The  nodes  also  need  to 
be  spaced  closely  enough  that  the  grid  resolution  becomes  acceptable,  particularly  for 
unsteady  or  multi-dimensional  heat  transfer. 

The  objective  of  all  thennal  analysis  codes,  such  as  SIND  A,  is  to  solve  the 
general  heat  equation  (Equation  2.4).  The  two  commonly  used  numerical  solution 
approaches  of  heat  transfer  problems  are  the  finite  difference  method  (FDM)  and  the 
finite  element  method  (FEM). 

The  FDM  approximates  the  physical  objects  and  makes  some  assumptions  such  as 
the  nodes  are  isothermal  and  physical  properties  are  constant  within  a  node.  Each  finite 
difference  node  is  located  at  the  center  of  the  mass.  The  heat  equation  is  converted  into  a 
system  of  finite  difference  equations  from  an  FDM  mesh.  The  FDM  is  based  on  the 
Taylor  series  approximation.  The  numerical  error  associated  with  the  calculation  can  be 
estimated  through  knowledge  of  the  order  of  the  Taylor  series  expansions. 

The  FEM  is  based  on  discretizing  the  structure  or  system  into  numerous  elements. 
Even  though  each  finite  element  model  includes  hundreds  of  elements  similar  to  finite 
difference  models,  finite  element  models  depend  on  geometry  and  application.  Element 
nodes  are  typically  located  at  the  elements’  corners.  Parameter  values,  such  as 
temperatures,  are  known  or  calculated  at  element  nodes.  Variations  within  the  element 
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are  calculated  by  interpolation  functions.  Thus  the  properties,  such  as  temperatures,  can 
vary  across  the  element. 

The  FDM  is  extensively  used  for  building  spacecraft-system  thermal  models. 
Heat  transfer  problems  that  are  primarily  driven  by  radiation  can  be  easily  solved  using 
the  FDM.  On  the  other  hand,  the  FEM  is  generally  used  in  structural  analysis  and  it  is 
excellent  for  solving  thermal/stress  problems.  For  example,  using  the  FEM  in  thennal 
analysis  of  rocket  nozzles  is  common. 

Energy  absorbed  on  surfaces  of  spacecraft  and  radiation  interchange  factors 
between  surfaces  are  calculated  by  codes  that  use  either  the  gray  diffuse  assumption  or 
the  Monte  Carlo  method.  The  codes  using  gray  diffuse  assumption  assume  that  the 
emittance,  absorptance,  and  transmittance  of  a  surface  are  independent  of  wavelength  and 
direction,  and  all  energy  from  a  surface  is  emitted  and  reflected  diffusely.  On  the  other 
hand,  codes  using  Monte  Carlo  technique  are  not  limited  to  these  assumptions.  Although 
Monte  Carlo  method  can  handle  complex  geometries,  large  detailed  thermal  models  can 
greatly  increase  the  execution  time  of  the  Monte  Carlo  software  [Gilmore,  2002]. 

Thermal  Desktop 

Autodesk’s  thennal  desktop  (TD)  is  a  thermal  modeling  program  that  allows  a 
user  to  build,  analyze,  and  post-process  thennal  models  [C&R  Technologies,  2013],  TD 
does  thennal  calculations  using  abstract  network,  finite  difference,  and  finite  element 
modeling  methods.  TD  runs  as  an  AutoCAD®  application  and  CAD  techniques  for 
generating  geometry  can  be  used  for  generating  thermal  models.  Monte-Carlo  ray-tracing 
algorithm  is  used  by  TD  to  compute  radiation  exchange  factors  and  view  factors. 
SINDA/FLUINT,  SINDA/G  and  ESATAN  are  thennal  analyzer  options. 
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SINDA/FLUINT  is  used  as  the  analyzer  in  this  research.  The  user  can  assign  optical 
properties  to  thermal  model  surfaces.  TD  also  provides  a  complete  external  heating 
environment  definition  and  viewing  facility,  which  allows  users  to  do  on-orbit  analyses. 
Detailed  usage  infonnation  about  TD  is  presented  in  Chapter  3. 

Spacecraft  Thermal  Testing 

In  a  spacecraft  thermal  management  process,  selected  thennal  design  is  validated 
and  thermally  successful  operation  of  the  spacecraft  is  ensured  by  ground  thermal  tests. 
These  tests  simulate  the  expected  flight  environments  for  the  spacecraft.  Test  facilities 
with  vacuum  are  used  to  minimize  non-flight  thermal  influences  and  simulate  on-orbit 
environments.  Test  categories  and  test  types  for  spacecraft  are  introduced  in  this  chapter 
based  on  Spacecraft  Thermal  Control  Handbook  by  Gilmore  [Gilmore,  2002]. 

There  are  three  test  categories:  development,  qualification,  and  acceptance  tests. 
Development  tests  are  used  to  validate  the  design  concepts  and  to  reduce  of  risk  in 
committing  designs  to  hardware  fabrication.  Development  test  requirements  are  unique  to 
test  objectives  and  are  not  specified  in  military  and  commercial  standards.  These  tests  are 
also  known  as  engineering  tests.  Qualification  tests  are  conducted  to  demonstrate  that  the 
design  and  manufacturing  process  meet  the  requirements.  Qualification  tests  simulate 
more  severe  conditions  than  expected  and  because  of  the  severity  of  the  test  enviromnent, 
qualification  hardware  is  typically  not  flown.  The  test  item  for  these  tests  is  produced 
from  the  same  drawings  using  the  same  materials.  Acceptance  tests  provide  quality 
control  assurance  against  workmanship  and  material  deficiencies.  These  tests  are 
conducted  after  qualification  tests  and  they  prove  the  flightworthiness  of  the  space 
system. 
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The  test  categories  explained  above  include  different  types  of  tests  such  as 
thennal  cycle  tests,  thermal  vacuum  tests,  burn-in  tests,  and  thennal  balance  tests. 
Thermal  cycle  tests  subject  the  test  article  to  a  number  of  cycles  of  hot  and  cold 
temperatures  in  an  ambient  air  or  gaseous  nitrogen  environment  whereas  thermal  vacuum 
tests  do  the  thermal  cycle  tests  in  a  vacuum  environment.  Thennal  cycle  tests’  primary 
purpose  is  to  screen  environmental  stress  and  reveal  workmanship  or  material  defects.  On 
the  other  hand,  thermal  vacuum  tests’  primary  purpose  is  performance  verification, 
because  these  are  conducted  without  convection  and  the  most  realistic  ground  simulation 
of  the  space  environment.  Burn-in  tests  are  typically  part  of  unit  thennal  cycle  tests  in 
which  additional  test  time  is  applied.  The  unit  is  either  held  or  cycled  at  an  elevated  hot 
temperature.  Thermal  balance  tests  have  two  purposes:  verification  of  the  thennal  control 
subsystem  and  correlation  of  thermal  analytic  models.  Correlation  within  1  and  5  °C  is 
typically  desired  depending  on  sensitivity  [Birur  et  ah,  2014].  Test  phases  simulate  hot 
and  cold  flight  conditions  to  gather  steady-state  temperature  data.  These  data  are 
compared  to  model  predictions  and  they  verify  all  aspects  of  the  thennal  hardware  and 
software. 

There  are  a  series  of  documents  which  specify  and  describe  the  requirements  for 
spacecraft  thermal  testing.  MIL-STD-1 540D  is  the  document  for  US  military  applications. 
It  was  published  in  1999.  This  document  provides  a  guideline  including  expected 
methodologies  and  acceptance  testing  requirements  without  specifically  directing  test 
practices  and  procedures  [Gilmore,  2002],  GSFC-STD-7000A,  General  Environmental 
Verification  Standard  (GEVS),  document  provides  requirements  and  guidelines  for 
environmental  verification  programs  and  describes  methods  for  implementing  those 
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requirements  for  NASA  applications  [NASA  GEVS,  2013].  GSFC-STD-7000A  directs 
thennal  test  practices  and  procedures  in  detail. 

Even  though  CubeSats  do  not  have  long  mission  lives,  thermal  tests  can  certainly 
increase  operational  lifetime  and  reduce  failure  rates  [Swartwout,  2013].  Various  tests, 
such  as  qualification  and  acceptance  tests,  are  typically  conducted  to  verify  the  thennal 
design  and  to  qualify  the  CubeSat  for  the  mission. 

CubeSat  Thermal  Design  Research  Overview 

Even  though  the  interest  in  CubeSats  has  increased  over  the  past  decade,  there  is 
little  publicly  released  research  about  thermal  design  and  analysis  of  the  CubeSats.  The 
following  paragraphs  present  some  of  the  researches  conducted  in  this  area. 

In  2004,  Czernik  at  the  University  of  Applied  Sciences  in  Aachen,  Germany,  did  a 
thennal  design  and  the  computational  analysis  of  a  1U  CubeSat  called  Compass- 1,  which 
is  designed  for  a  sun-synchronous  LEO  orbit  [Czemik,  2004].  Active  thermal  control 
hardware,  such  as  heaters,  was  included  as  well  as  passive  hardware  in  the  design. 
ANSYS  software,  which  uses  FEM  to  calculate  nodal  temperatures,  was  chosen  for 
computational  analysis,  but  experimental  analysis  was  not  performed  to  validate 
computational  analysis  in  the  research. 

In  2009,  Diaz-Aguado  et  al.  at  NASA  Ames  Research  Center  in  Moffett  Field, 

California  designed  the  thermal  control  system  of  3U  CubeSat  named  PharmaSat 

[Diaz-Aguado  et  al.,  2009].  The  thennal  design  included  active  and  passive  thermal 

control  hardware  such  as  heaters,  MLI,  and  low  thermal  conductance  materials. 

Computational  analysis  was  perfonned  in  TD  and  experimental  tests  conducted  in  a 

thermal  vacuum  chamber.  Worst  hot  case  and  worst  cold  case  results  of  the  experimental 
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tests  were  validated  with  computational  simulations  with  +5  °C  error  acceptance. 
However,  Diaz-Aguado  et  al.  provide  little  information  about  the  thermal  model. 

Elhady,  at  National  Authority  for  Remote  Sensing  &  Space  Science  NARSS  in 
Cairo,  Egypt,  did  a  thermal  design  and  analysis  of  a  740  mm  x  740  mm  x  722  mm  micro 
satellite  in  2010  [Elhady,  2010],  The  effect  of  thennal  contact  resistance  on  the  satellite 
was  also  analyzed.  An  in-house  developed  software  package  was  used  for  computational 
analysis.  However,  the  research  did  not  include  the  experimental  analysis. 

Smith  conducted  a  thennal  analysis  of  a  CubeSat  called  NPS-SCAT  at  the  Naval 
Postgraduate  School  (NPS)  in  2011  [Smith,  2011],  For  computational  analysis,  NX-6 
I-DEAS  software  was  used  with  FEM  for  predictions.  Experimental  test  results  of  TVAC 
were  used  to  validate  computational  models.  There  were  significant  temperature 
differences  during  the  cold  soak  caused  by  some  effects  such  as  simplification  of  the 
model  compared  to  the  complexity  of  the  satellite  and  inaccuracies  in  the  TVAC 
temperature  profile.  These  effects  were  discussed  in  research  and  some  recommendations 
were  provided  for  future  work. 

In  2012,  thermal  design  and  computational  analysis  of  a  1U  CubeSat  were 
performed  by  Dinh  at  San  Jose  State  University  in  California  [Dinh,  2012].  Dinh  used  TD 
and  ANSYS  Icepak  as  thermal  software;  however,  the  research  did  not  include 
experiments  for  validation  of  the  models. 

Different  from  all  researches  discussed  above,  this  presented  thesis  provides  a 
modular  thermal  design  for  the  CubeSats.  This  design  was  created  using  both 
experimental  and  computational  methodologies.  All  computational  models  were 
validated  using  experimental  data  in  this  thesis  effort. 
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Finally,  Hatzung’s  research  is  briefly  discussed  in  this  chapter  to  provide  a 
background  for  the  first  task,  provided  on  page  3,  of  this  presented  research  [Hatzung, 
2014].  At  AFIT  in  2013,  Hatzung  used  TD  to  create  a  computational  thermal  model  of 
the  AFIT’s  TVAC.  The  purpose  of  the  TVAC  model  is  to  validate  AFIT’s  future  CubeSat 
thennal  models  by  comparing  with  experimental  test  results.  Hatzung  did  experiments  to 
validate  his  model  using  a  10  in  x  10  in  aluminum  plate  and  demonstrated  that  the  model 
is  valid  for  use  in  steady  state  temperature  predictions  for  model  validation.  Due  to  time 
considerations,  he  could  not  conduct  experiments  with  another  test  article  for  greater 
validation. 

Summary 

This  chapter  provided  a  theoretical  basis  through  an  overview  of  heat  transfer  and 
Earth  orbit  thennal  environment.  Thennal  design  process  and  commonly  used  thermal 
control  hardware  were  described.  As  emphasized  before,  the  design  process  is  a 
combination  of  design  selection  and  validation  tests.  Both  experimental  and 
computational  thennal  analyses  are  crucial  parts  of  the  thermal  design  process. 
Computational  and  experimental  tests  provide  validation  not  only  for  the  selected  design 
but  also  for  each  other.  Finally,  to  contribute  to  basic  understanding  of  the  thennal  design 
and  analysis  effort  of  CubeSats  in  the  past  decade,  a  research  overview  was  provided. 
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ITT.  Methodology 


The  purpose  of  this  chapter  is  to  describe  the  methodologies  used  in  this  research. 
First,  the  experimental  test  system  of  the  thesis,  Abbess  produced  TVAC  test  system,  is 
introduced.  Second,  an  overview  of  Thermal  Desktop  software  is  provided  as  the 
computational  methodology  and  a  brief  thennal  model  design  and  test  process  is 
presented  to  provide  a  basic  understanding  of  the  software  logic.  Third,  test  articles  and 
procedures  are  discussed.  Fourth,  validation  process  is  discussed.  Finally,  the  heat 
dissipater  design  process  is  presented. 

Experimental  Test  Methodology 

AFIT’s  TVAC  was  used  for  experimental  tests  in  this  research.  The  TVAC 
system  is  described  and  a  brief  operation  process  is  introduced  below. 

TVAC 

TVAC  approximates  specific  thermal,  vacuum  pressure,  and  illumination 
conditions  in  order  to  observe  articles  under  test,  which  may  experience  similar 
conditions  in  actual  use  [Thermal/Vacuum  Solar,  2012],  Vacuum  and  temperature 
conditions  are  established  using  on  board  controls.  These  controls  activate  thermal 
system  components  and  turbo  molecular  pumping  system  to  set  the  desired  conditions 
within  the  chamber.  A  front  view  of  the  TVAC  is  shown  in  Figure  2.  Some  of  the  major 
components  of  the  TVAC  are  introduced  below. 
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Figure  2:  Front  View  of  TVAC 

The  vacuum  chamber  is  a  front-loading,  insulated  stainless  steel  chamber,  which 
serves  to  condition  the  chamber  for  high  vacuum  performance  down  to  1x10°  Torr. 
Temperature  within  the  chamber  is  measured  by  built-in  temperature  sensors,  which  are 
located  on  the  platen  and  the  shroud.  The  chamber  has  process  ports,  which  were  used  for 
thermocouple  and  power  connections  in  experiments. 

The  thermal  platen/shroud  assembly  consists  of  a  copper  thermal  platen  and  a 
shroud  mounted  in  vacuum  chamber  with  engaged  heating/cooling  loop  as  shown  in 
Figure  3.  The  platen  and  shroud  facilitate  heating/cooling  profiles  (— 40°C  to  125°C  )  in 
the  chamber  under  vacuum.  These  profiles  are  controlled  by  an  integrated  control 
interface  provided  by  a  touch  screen  PC.  The  platen  is  gold  plated  and  the  shroud  is 
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Aeroglaze  black  paint  coated.  They  are  each  heated  and  cooled  individually  by  an 
engaged  thennal  media  loop  piping  from  a  thermal  bath  circulator  system  manufactured 
by  Thermo  Fisher  Scientific.  The  front  wall  of  the  shroud  is  removable  and  is  not  actively 
heat  controlled.  Conduction  with  the  other  walls  of  the  shroud  mainly  determines  the 


temperature  of  the  front  wall  during  test  profiles. 


Thennal  Loop  Piping 


Gold  Plated  Copper 
Thermal  Platen 


Stainless  Steel 
Vacuum  Chamber 


Vacuum  Gauges 


Opening  to  Solar 
Simulator  Gate 


Aeroglaze  Coated 
Copper  Shroud 


Figure  3:  Open  Door  View  of  Chamber 

The  vacuum  system  consists  of  a  roughing  pump  and  a  turbo  pump.  The  roughing 
pump  removes  the  bulk  of  gas  from  the  chamber  from  atmosphere  down  to  a  mid-level 
vacuum  range.  The  turbo  pump  is  a  6,000  to  49,000  RPM  pump  used  to  achieve  high  and 
ultra  high  vacuum  pressure  in  chamber  (P  <  lxlO"5  Torr).  This  pump  system  is  powered 
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by  a  digital  control  unit  with  pressure  data  provided  by  two  separate  digital  vacuum 


gauges. 


The  solar  simulator  is  designed  to  achieve  one  sun  equivalent  illumination  with 
three-degree  collimation  over  a  12  in  x  12  in  area.  A  1.6  kW  lamp  is  used  for  the 
illumination  and  managed  by  an  arc  lamp  power  supply  and  digital  exposure  controllers. 
The  solar  simulator  was  unnecessary  for  this  research. 


Figure  4:  Front  View  of  Thermal  Vacuum  Control  Enclosure  and  Main  Screen 

The  TVAC  is  controlled  by  an  integrated  control  interface  provided  by  a  touch 
screen  PC  as  expressed  above.  The  main  screen  of  PC  is  as  shown  in  Figure  4.  Real  time 
plate,  shroud,  and  bake-out  control  and  limit  temperatures  are  displayed  on  this  screen. 
On  the  left,  chamber  pressure,  in  Torr,  is  displayed  according  to  two  digital  vacuum 
gauges,  which  are  shown  in  Figure  4.  A  user  is  able  to  navigate  the  Pressure  Control, 
Temperature  Control,  and  Solar  Simulator  interfaces  as  well  as  the  Setup  and  Relay 
Status  screens  using  the  buttons  on  main  screen.  The  user  can  also  create  or  select  a  file 
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to  record  experiment  data  using  buttons  at  the  top  center  of  the  main  screen.  The  data  file 
is  a  text  file,  which  can  be  opened  in  Notepad  or  Excel  format. 

To  conduct  an  experiment,  when  the  chamber  is  loaded  with  a  properly  prepared 
test  article,  the  user  secures  all  chamber  doors  and  decreases  the  pressure  using  the 
Start/Stop  button  on  the  pressure  control  screen  shown  in  Figure  5.  On  the  left,  chamber 
pressure  in  Torr  is  displayed  according  to  two  digital  vacuum  gauges  in  digital  data 
format.  On  the  right  half  of  the  screen,  a  graphic  display  charting  chamber  pressure 
against  time  is  shown  with  color-coded  traces  of  the  two  gauges  in  real  time. 


Figure  5:  Pressure  Control  Screen 

When  the  chamber  is  under  vacuum,  the  user  may  individually  enter  a 
temperature  value  for  plate  and  shroud  using  the  temperature  control  screen,  which  is 
shown  in  Figure  6.  The  center  of  this  screen  is  devoted  to  graphic  display  charting 
temperature  vs.  time  with  color-coded  trace  of  plate,  shroud,  and  bakeout  temperatures  in 
real  time. 
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Figure  6:  Temperature  Control  Screen 

When  the  experiment  is  finished,  the  user  can  start  the  controlled  depressurization 
process  using  the  Start/Stop  button  on  pressure  control  screen  after  returning  the  chamber 
to  room  temperature. 

Computational  Simulation  Methodology 

As  discussed  in  the  previous  chapter,  Thermal  Desktop  (TD)  was  used  for 
computational  simulations  in  this  research.  Even  though  the  manual  [C&R  Technologies, 
2013],  which  is  available  online,  provides  very  detailed  information,  a  brief  thermal 
model  design  and  test  process  is  presented  here  to  assist  reader  in  future  discussion. 

Thermal  Desktop 

C&R  Thermal  Desktop®  is  a  program  that  allows  the  user  to  build,  analyze,  and 
post-process  thermal  models  [C&R  Technologies,  2013].  TD  runs  as  an  AutoCAD 
application  and  can  import  some  other  models,  such  as  TRASYS,  NEVADA,  TSS,  and 
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NASTRAN,  and  I-DEAS.  Every  new  object  is  associated  with  a  layer,  and  layers  are 
used  to  organize  the  design  as  in  AutoCAD. 

TD  objects  are  generated  using  thermal-specific  surface  types  such  as  rectangle, 
disc,  cylinder  and  solids  such  as  brick,  solid  cylinder.  CAD  surfaces  created  by  built-in 
CAD  construction  techniques  can  be  converted  to  TD  surfaces  as  well.  Starting  the 
drawing  process  with  defining  materials  and  their  thermophysical  and  optical  properties 
may  be  timesaving  for  a  modeling  process.  Units  of  thermophysical  properties  can  be 
changed  by  the  user.  The  SI  unit  system  (W/m/K  for  conductivity,  J/kg/K  for  specific 
heat,  and  kg/m3  for  density)  was  used  for  this  thesis  effort. 

When  a  surface  or  solid  is  created,  the  Thin  Shell  Data  dialog  box  appears  as 
shown  in  Figure  7. 


Figure  7:  TD  Thin  Shell  Data  Dialog  Box 

The  user  selects  the  node  type  and  the  number  of  the  nodes  in  each  direction  from  the 
Subdivision  tab.  The  user  selects  Centered  Nodes  button  for  finite  difference  objects  and 
Edge  Nodes  button  for  finite  element  objects.  If  the  processing  speed  and  minimizing  the 
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number  of  the  nodes  are  important  (for  large  models),  the  user  may  consider  finite 
difference  objects  when  the  shape  of  the  geometry  fits  one  of  the  built-in  finite  difference 
surfaces  and  solids.  When  creating  complex  models  with  irregular  curves,  finite  element 
objects  are  preferable.  However,  the  user  should  expect  increased  processing  time  with 
increased  complexity  of  the  model  when  using  finite  element  objects.  Furthermore,  care 
must  be  taken  when  selecting  edge  or  centered  nodes,  particularly  in  the  case  of  creating 
conduction  between  separate  surfaces  using  a  contactor,  which  is  described  below.  Yhe 
contactors  do  not  calculate  the  conduction  within  any  nodal  breakdown  in  an  object. 
When  conduction  created  from  the  lower  two  objects,  shown  in  Figure  8,  to  upper  object, 
the  contactor  functionality  calculates  the  conduction  from  the  lower  objects’  surfaces  to 
upper  object’s  near  surface.  This  means,  in  Figure  8,  conductions  between  Point  1  and 
Point  2  (conduction  within  any  nodal  breakdown)  are  not  calculated.  If  that  conduction 
is  important,  edge  nodes  should  be  used.  Using  small  finite  element  objects  in  a  large 
finite  difference  model  does  not  create  a  serious  effect  on  the  precision  of  the  model.  A 
few  objects  were  created  with  edge  nodes,  such  as  circuit  board  and  1U  stack  rods,  in  this 
research. 
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Figure  8:  Contactor  Example 
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The  object  material  is  selected  using  the  Cond/Cap  tab.  Optical  properties  of  the 
surfaces  are  selected  on  the  Radiation  tab  for  radiation  calculations.  All  surfaces  and 
solids  can  be  modified  (moving,  rotating,  stretching  etc.)  by  using  built-in  AutoCAD 
commands. 

Since  no  fluid  was  used  for  heat  dissipation,  convective  heat  transfer  was  ignored 
in  this  research.  TD  uses  optical  properties  and  physical  locations  of  the  objects  to 
calculate  radiation.  Nevertheless,  conductive  heat  transfer  should  be  created  by  the  user 
on  TD  using  conductors,  contactors  or  merge  nodes. 

A  conductor  can  be  created  using  one  of  the  three  choices:  Node-to-Node,  Node- 
to-Nodes,  or  Node-to-Surface  Conductor.  In  this  conduction  type,  the  user  should  select  a 
single  node  to  detennine  the  from  object  of  conduction.  Conductors  allow  specifying  the 
conductivity  value  of  the  intervening  material  as  a  function  of  time  or  temperature 
difference.  Since  heat  transfer  calculations  from  an  individual  node  were  not  necessary  in 
this  research,  the  conductor  was  not  use  as  a  conduction  type. 

The  contactors  are  created  to  model  conduction  between  two  or  more  high  level 
objects,  such  as  2D  surfaces  or  3D  solids.  For  finite  conductions  between  surfaces  and 
solids,  the  contactors  are  the  best  option  and  widely  used  in  this  thesis.  The  contactor 
divides  the  conductive  heat  transfer  into  a  k-A/L  or  tv  A  tenn,  where  k  is  the  conductivity, 
A  is  the  area,  L  is  the  thickness  of  the  interface  and  h  is  the  heat  transfer  coefficient.  To 
create  the  contactor,  the  user  can  provide  a  fixed  heat  transfer  coefficient  and  the 
contacted  area  is  calculated  by  the  contactor.  Alternatively,  use  of  the  material  option 
allows  the  user  to  choose  an  intervening  material.  With  this  option,  the  contactor  takes 
the  conductivity  value  from  the  thermophysical  properties  and  calculates  the  area  values 
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from  the  geometry.  The  user  needs  only  to  input  the  thickness  of  the  interface.  The 
contactor  also  allows  the  user  to  pick  the  Ray  Trace  Algorithm,  which  shoots  rays 
perpendicular  to  the  “from  surface”  and  deems  the  closest  node  to  be  whatever  node  the 
ray  first  encounters.  Since  this  is  a  very  fast  algorithm,  it  was  used  in  this  research  for  the 
contactors  created  between  parallel  surfaces. 

For  very  large  conductances  value  or  perfect  contact,  merging  nodes  should  be 
considered.  Merging  nodes,  where  possible,  increases  the  processing  speed  of  the  model 
instead  of  creating  the  contactors  with  large  conductances.  As  an  example,  nodes  at  the 
edges  are  merged  when  the  aluminum  chassis  of  the  CubeSat  was  created  using  a  small 
piece  of  aluminum  solids  to  simulate  large  conductances.  However,  the  contactors  were 
used  to  create  conduction  between  the  chassis  and  its  covers  to  simulate  relatively  smaller 
conductances. 

When  the  model  is  created,  using  the  techniques  explained  above,  the  user  can 
simulate  an  orbit  using  Heating  Rate  Case  Manager  dialog  box  as  shown  in  Figure  9. 


Heating  Rate  Case  Manager 


Current  Heating  Rate  Case: 


[  Display  Orbit  ]  |  Done  | 


_ 0 


|  Add  |  |  Import  | 


|  Help  ] 


Figure  9:  Heating  Rate  Case  Manager  Dialog  Box  and  Orbit  Dialog  Box 
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The  Create  New  External  Heating  Environment  dialog  box  appears  when  selecting  the 
Add  button.  TD  allows  the  user  to  define  many  external  heating  environments  such  as 
basic  orbit,  Keplerian  orbit,  and  free  molecular  heating  with  a  reference  orbit.  A  basic 
orbit  was  used  in  this  research  to  get  the  orbit  performance  of  the  heat  dissipater  design. 
The  basic  orbit  constellation  allows  the  user  to  define  the  altitude  and  beta  angle  of  a 
circular  orbit  as  shown  in  Figure  9.  The  Orbit  dialog  box  also  allows  user  to  specify  the 
orientation  of  the  spacecraft,  the  locations  (positions)  in  orbit  for  which  heating  rates  will 
be  calculated,  planetary  data  which  the  satellite  orbits  such  as  earth,  sun,  moon,  and 
orbital  heating  rate  data  such  as  direct  solar,  albedo  and  IR  emittance  of  the  planet.  The 
orbit  and  selected  positions  can  be  displayed  as  shown  in  Figure  10. 


Figure  10:  Orbit  Display 

Color-coded  coordinate  systems  (XYZ  axes)  on  the  orbit  display  depict  the  spacecraft’s 


orientation  at  each  position  with  respect  to  the  celestial  coordinate  system,  which  is  at  the 
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center  of  the  Earth  (X-red,  Y-green,  and  Z-blue).  Purple  cylinders  show  the  shadow  zone 
of  the  Earth. 

After  defining  the  orbit  and  the  satellite  orientation,  the  user  can  compute  heating 
rates  for  solar  albedo  and  IR  planetshine  using  a  Monte  Carlo  method.  The  user  can  also 
compute  orbital  heating  rates  individually  or  get  the  error  rate  in  calculations  by  using  Set 
HR  Dataset  Properties  dialog  box  as  a  post-processing  as  shown  in  Figure  11. 


Figure  11:  Set  HR  Dataset  Properties  Dialog  Box 

When  the  model  and  external  heating  environment  (orbit)  are  created,  the  Case 
Set  Manager  is  used  to  conduct  a  simulation.  The  Case  Set  Manager’s  primary  purpose  is 
to  set  up  different  thermal  analysis  cases  and  to  make  temperature  calculations  from 
radiation  calculations.  First,  it  performs  the  radiation  and  heating  rate  tasks  of  the  current 
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case,  computes  and  outputs  the  node  temperature  values  and  node -to  node  radiation 
conductor  data,  builds  and  runs  a  SINDA  model,  and  finally  displays  the  temperature 
results  mapped  on  the  thermal  model.  A  case  can  be  added  or  the  current  case  can  be 
edited  by  Case  Set  Manager  dialog  box.  The  Editing  Case  Set  dialog  box  appear  as 
shown  in  Figure  12,  when  add  or  edit  button  is  selected  on  Case  Set  Manager  dialog  box. 


Figure  12:  Editing  Case  Set  Dialog  Box 

New  radiation  tasks  can  be  added  or  the  current  ones  can  be  edited  on  the  Radiation 
Tasks  tab.  The  user  defines  solution  types  such  as  transient  or  steady  state  on  the 
Calculations  tab  and  selects  type  of  output  data  on  Output  tab.  The  selected  type  of  data 
is  written  on  the  output  files  such  as  temperatures,  capacitances  or  radiation  data  (e.g., 
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radiation  calculation  results  such  as  number  of  rays  shot  for  each  node  and  weighted 
error). 

When  the  calculations  are  performed,  calculation  results  can  be  examined  by  post¬ 
processing  such  as  color  post-processing  (displaying  calculation  values  on  graphical 
object  as  colors),  X-Y  plotting  (as  a  function  of  time),  query  node  (by  selecting  graphical 
objects),  and  results  queries  (generating  summary  fdes  based  on  one  or  more  result  fdes). 
Temperature  color  post-processing  figures  are  provided  for  the  computational  simulations 
conducted  in  this  research  in  Chapter  4. 

Test  Articles  and  Experiments 

A  simple  test  article,  a  6U  CubeSat  Chassis,  was  used  for  the  initial  tests  to 
achieve  the  first  task  of  the  thesis,  which  is  to  provide  further  validation  for  the 
previously  generated  TVAC  environment  computational  model  [Hatzung,  2014].  The 
Chassis  is  entirely  made  from  T-6061  aluminum.  The  Chassis  was  instrumented  with 
thennocouples  and  hung  from  the  test  stand  using  a  monofilament  fishing  line  as  shown 
in  Figure  13  to  minimize  conductive  heat  transfer. 


38 


(a)  Side  View 


Figure  13:  Two  Views  of  the  Chassis  Hung  from  the  Test  Stand 

Two  circuit  boards  with  resistors  were  used  in  experiments  in  this  research.  The  first  one 
was  a  custom  PC  104  circuit  board  with  ten  resistors,  as  shown  in  Figure  14  (a),  which 
was  named  as  ten-heater  circuit  board  (THCB).  Ten  Vishay  WSC4527  wirewound 
resistors,  each  capable  of  dissipating  2W  of  power,  were  mounted  on  the  Isola  FR402 
circuit  board  to  simulate  processors  in  a  CubeSat.  The  second  circuit  board  was  built  with 
a  Vishay  10W  Resistor  as  shown  in  Figure  14  (b)  and  named  as  single-heater  circuit 
board  (SFICB)  for  this  research. 
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(a)  The  Ten-Heater  Circuit 
Board, THCB 


(b)  The  Single-HeaterCircuit 
Board,  SHCB 


Figure  14:  Circuit  Boards  with  Resistors 

The  chassis  was  covered  with  aluminum  plates  for  succeeding  experiments  to 
simulate  a  fully  enclosed  6U  CubeSat  as  shown  in  Figure  15.  The  covered  chassis 
including  heat  sources  was  named  as  sample  CubeSat  in  this  research.  The  covers  were 
made  of  T-6061  aluminum.  A  small  gap  was  cut  on  the  cover  to  allow  thermocouple  and 
power  connections  between  the  TVAC  and  CubeSat  interior  as  shown  in  bottom  right 
comer  of  Figure  15. 


Figure  15:  The  Sample  6U  CubeSat 
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All  experiments  are  named  with  a  letter  to  reduce  confusion  in  this  research.  For 
example,  chassis-only  experiment,  which  is  discussed  below,  is  named  as  experimental 
test- A  (ET-A).  In  case  of  an  experiment  is  also  modeled  computationally  in  TD,  the 
computational  simulation  is  also  named  with  the  same  method.  For  example,  chassis-only 
computational  simulation  is  named  as  computational  simulation-A  (CS-A).  All 
experiments  and  computational  simulations  conducted  in  this  research  are  presented  in 
Table  1.  Even  though  many  items  in  Table  1,  such  as  radiators  or  heat  dissipater 
configuration,  have  not  been  introduced  at  this  point,  the  table  is  presented  before 
experiment  and  simulation  discussions  to  give  the  reader  a  big  picture. 
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Table  1:  Test  Specifications  and  Explanations 


Experimental  Test 
Name 

Computational 
Simulation  Name 

Explanation 

Purpose 

Heater 

Heat 

Dissipater 

Radiator  Size 

TVAC  Test  Point 
(°C) 

Power  (W) 

ET-A 

CS-A 

Chassis-Only 

Provide  greater  validation 
for  previously  created 
TVAC  computational 
model 

- 

- 

- 

-10 

30 

60 

- 

ET-B 

Chassis  with 

Heaters 

Acquire  thermal  behavior 
of  the  heaters  in  vacuum 

THCB 

- 

- 

-15 

1.8 

2.7 

4.5 

6.3 

ET-C 

Sample 

CubeSat 

Power  Level 

Determination 

Power  level  determination 
for  subsequent  experiments 

THCB 

- 

- 

-15 

1.8 

2.7 

4.5 

ET-D 

CS-D 

Sample 

CubeSat 

Acquire  thermal  behavior 
of  the  sample  CubeSat  and 
computationally  model  it 

THCB 

- 

- 

-15 

1.8 

4.5 

ET-E 

Sample 

CubeSat 

Determine  thermal  behavior 

of  test  article  with  heat 
dissipater  box  and  radiator 
when  there  is  no  connection 

between  the  box  and  the 

radiator 

THCB 

Box  and 
Radiator  (No 
thermal  strap 
connection) 

1U  (Bare  Copper) 

-15 

1.8 

4.5 

ET-F 

Sample 

CubeSat 

Determine  the  performance 
of  the  heat  dissipater  with 
bare  copper  radiator  (no 
surface  finish) 

THCB 

Box,  straps 
and  radiator 
(Poor 

Conduction 

in  Braid 
Connections) 

1U  (Bare  Copper) 

-15 

1.8 

4.5 
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Experimental  Test 
Name 

Computational 
Simulation  Name 

Explanation 

Purpose 

Heater 

Heat 

Dissipater 

Radiator  Size 

TVAC  Test  Point 
(°C) 

Power  (W) 

ET-G 

Sample 

CubeSat 

Determine  the  effect  of 

surface  finish  of  the 

radiator  in  the  heat 
dissipater  design  (with 
comparison  of  ET-F) 

THCB 

Box,  straps 
and  radiator 
(Poor 

Conduction 

in  Braid 
Connections) 

1U  (Black  Tape) 

-15 

1.8 

4.5 

ET-H 

CS-H 

Sample 

CubeSat, 

Basic 

Configuration 
of  the  Heat 
Dissipater 

Determine  the  importance 
of  thermal  strap  connection 
quality  in  the  heat  dissipater 
design  (with  comparison  of 
ET-G) 

THCB 

Basic 

Configuration 
(Thermal 
Straps  with 
High 

Temperature 

Soldering) 

1U  (Black  Tape) 

-15 

1.8 

4.5 

11 

ET-I 

Sample 

CubeSat 

Determine  the  effect  of 

radiator  size  in  the  heat 
dissipater  design 

THCB 

Basic 

Configuration 

3U  (Bare 
Copper) 

-15 

1.8 

4.5 

ET-J 

Sample 

CubeSat 

Determine  the  effect  of 

radiator  size  and  radiator 

surface  finish  in  the  heat 
dissipater  design 

THCB 

Basic 

Configuration 

3U  (Black 
Tape) 

-15 

1.8 

4.5 

7.2 

ET-K 

Sample 

CubeSat 

Determine  the  effect  of  a 
deployable  radiator  in  the 
heat  dissipater  design 

THCB 

Basic 

Configuration 

1U  Side  + 
Deployable 

-15 

1.8 

4.5 

ET-L 

Sample 

CubeSat 

Acquire  the  thermal 
behavior  of  the  CubeSat 

with  SHCB 

SHCB 

- 

- 

-15 

3.3 
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Experimental  Test 
Name 

Computational 
Simulation  Name 

Explanation 

Purpose 

Heater 

Heat 

Dissipater 

Radiator  Size 

TVAC  Test  Point 
(°C) 

Power  (W) 

ET-M 

Sample 

CubeSat 

Determine  the  effectiveness 
of  the  heat  dissipater 
design,  when  there  is 
indirect  contact  between  the 

heat  source  and  the  heat 
dissipater 

SHCB 

Full 

Configuration 
(Basic  + 
Aluminum 
Edges) 

1U  (Black  Tape) 

-15 

3.3 

ET-N 

Sample 

CubeSat 

Determine  the  effectiveness 
of  the  basic  configuration 
of  heat  dissipater,  in  case 
of  an  concentrated  heat 

source 

SHCB 

Basic 

Configuration 

1U  (Black 
Tape) 

-15 

3.3 

ET-0 

Sample 

CubeSat 

Determine  the 

effectiveness  of  the  heat 
dissipater’ s  box  shape  in 
comparison  with  ET-H 

THCB 

Basic 

Configuration 
w/o  Bottom 

Face 

1U  (Black 
Tape) 

-15 

1.8 

4.5 

CS-P 

0°  Beta  Angle 
Orbit 

Acquire  the  thermal 
behavior  of  the  sample 
CubeSat  in  the  cold  case 

orbit 

THCB 

- 

- 

- 

4.5 

CS-Q 

0°  Beta  Angle 
Orbit 

Determine  the  effectiveness 
of  the  heat  dissipater  in  the 
sample  cold  case  orbit 

THCB 

Basic 

Configuration 

1U  (Black 
Tape) 

- 

4.5 

13.5 

CS-R 

0°  Beta  Angle 
Orbit 

Determine  the  effect  of  the 
radiator  size  in  the  sample 
cold  case  orbit 

THCB 

Basic 

Configuration 

3U  (Black  Tape) 

- 

4.5 

19 
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Experimental  Test 
Name 

Computational 
Simulation  Name 

Explanation 

Purpose 

Heater 

Heat 

Dissipater 

Radiator  Size 

TVAC  Test  Point 
(°C) 

Power  (W) 

cs-s 

90°  Beta 
Angle  Orbit 

Acquire  the  thermal 
behavior  of  the  sample 
CubeSat  in  the  sample  hot 
case  orbit 

THCB 

- 

- 

4.5 

CS-T 

90°  Beta 
Angle  Orbit 

Determine  the  effectiveness 
of  the  heat  dissipater  in  the 
sample  hot  case  orbit 

THCB 

Basic 

Configuration 

1U  (Black 
Tape) 

- 

4.5 

15.5 

CS-U 

90°  Beta 
Angle  Orbit 

Determine  the  effect  of  the 
radiator  size  in  the  sample 
hot  case  orbit 

THCB 

Basic 

Configuration 

3U  (Black 
Tape) 

- 

4.5 

22 

In  experiments,  dT/dt  <  0.5  °C  /  hr  was  determined  as  thermal  equilibrium  and 
steady  state  criterion  as  explained  in  Chapter  1 .  Even  though  the  thermal  equilibrium  was 
achieved  approximately  in  15-18  hours  in  experiments,  all  experiments  were  conducted 
24  hours  to  provide  greater  confidence  to  data.  Due  to  some  failures  of  the 
thennocouples,  i.e.,  detachment  from  the  surface,  some  experiments  were  conducted 
repeatedly  to  acquire  dependable  data  throughout  this  research.  As  an  example,  in  some 
experiments  two  thennocouples  were  used  to  acquire  temperatures  from  the  circuit  board 
top  surface,  expecting  to  show  the  same  results,  because  the  circuit  board  has  roughly 
uniform  heat  dissipation  on  its  surface  (uniformity  of  the  heat  dissipation  of  the  circuit 
board  surface  will  be  discussed  in  the  next  chapter).  When  different  results  were  collected 
from  the  thermocouples,  tests  were  repeated.  Furthennore,  thennal  images  of  the  circuit 
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board  were  acquired  when  the  heaters  are  active  to  ensure  that  there  is  a  uniform  heat 
distribution  on  the  circuit  board.  Since  TVAC  does  not  allow  use  of  a  thennal  camera  due 
to  lack  of  an  observation  window,  IR  images  were  acquired  in  ambient  pressure  and 
temperature. 

Baseline  Experiments  and  Simulations 

Experiments  conducted  and  simulations  created  to  acquire  a  validated  CubeSat 
model  were  named  as  baseline  in  this  research.  Baseline  experiment  and  simulations 
consist  of  chassis-only  experiment  (ET-A),  chassis-only  simulation  (CS-A),  chassis  with 
heaters  experiment  (ET-B),  CubeSat  with  heater  experiment  (ET-C),  sample  CubeSat 
experiments  (ET-D),  and  sample  CubeSat  simulations  (CS-D).  Even  though  the  test 
article  is  the  same  in  ET-B  and  ET-C,  they  were  named  separately,  because  the  purposes 
of  the  experiments  were  different  from  each  other. 

Chassis-Only  Experiment  and  Simulation  (ET-A  and  CS-A) 

The  purpose  of  the  chassis  only  experiment  and  simulation  was  to  provide  further 
validation  for  the  previously  generated  TVAC  environment  computational  model  using  a 
6U  CubeSat  chassis  [Hatzung,  2014].  The  chassis  only  experiment  (ET-A)  was 
conducted  with  three  steady  state  temperature  set  points,  -10,  30,  and  60  °C.  Set  points 
were  selected  within  the  TVAC’s  temperature  envelope  to  be  able  to  perform  the  TVAC 
model  validation  at  cold,  hot,  and  room  temperatures.  Five  thennocouples  were  used  to 
obtain  temperature  data  from  the  chassis. 

In  the  chassis  only  computational  simulation  (CS-A),  an  analogous  computational 
model  of  the  chassis  was  created  in  TD  and  located  in  TVAC  computational  model  at  the 
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same  position  and  orientation  with  ET-A.  Photographs  and  TD  screenshots  of  ET-A  and 
CS-A  are  shown  in  Figure  16. 


(a)  The  Chassis,  Hung 
from  the  Test  Stand 


(b)  Thermal  Model  of  the 
Chassis,  Hung  from  the  Test 
Stand 


(c)  The  Chassis  inside 
TVAC 


(d)  Thermal  Model  of 
the  Chassis  inside  TVAC 


Figure  16:  Experiment  Pictures  and  TD  Screenshots  of  ET-A  and  CS-A 

Chassis  with  Heaters  Experiment  (ET-B)  and  Sample  CubeSat  Power  Level 
Determination  Experiment  (ET-C) 

The  purpose  of  ET-B  was  to  determine  the  thermal  behavior  of  the  heaters  under 
vacuum.  The  purpose  of  ET-C  was  to  determine  two  power  levels  for  subsequent 
experiments  to  reduce  the  experiment  time.  Additionally,  the  data  collected  from  ET-C 
was  compared  with  ET-B  and  this  information  was  used  to  examine  the  thermal  behavior 


47 


of  the  heaters  when  their  open  view  to  TV  AC  shroud  was  obstructed  with  aluminum 
plates. 

The  circuit  board  with  the  heaters  was  connected  to  an  external  power  supply 
through  a  feed-through  port  of  TV  AC.  The  power  level  was  set  using  the  external  power 
supply.  Since  the  cables  between  the  power  supply  and  the  circuit  board  were  relatively 
long  and  had  electrical  resistance,  power  loss  in  the  cables  was  measured  and  calculated 
as  10%.  For  this  reason,  the  power  levels  dissipated  at  the  circuit  board  were  used  to 
name  the  tests,  even  though  the  power  levels  measured  at  the  power  source  were  higher 
(i.e.,  the  power  at  the  power  source  was  2W  in  1.8W  tests;  and  it  was  5W  in  4.5W  tests). 
Thus,  the  THCB  model  was  run  with  the  power  matching  the  power  applied  on  the 
heaters,  not  the  power  supplied  by  the  power  supply.  Figure  17  shows  the  photographs 
ET-B. 


(a)  1U  Stack  including 
the  Circuit  Board 

(b)  The  Chassis  including 
Heaters  in  TVAC 

Figure  17:  Photographs  of  ET-B 

In  chassis  with  heaters  experiment,  the  TVAC  shroud  temperature  was  set  to  - 

15  °C  and  the  power  level  of  the  heaters  was  gradually  increased.  Experiments  were 

48 


started  with  1.8W  power  level,  and  the  power  was  increased  to  2.7W,  4.5W,  and  6.3W 
allowing  sufficient  time  for  steady  state  to  be  achieved  at  each  power  level.  The  top  face 
of  the  circuit  board  was  measured  at  each  test  point.  Experiment  terminated  at  the  6.3W 
power  level,  which  had  a  temperature  of  92.2  °C  on  the  circuit  board.  Since  the  maximum 
operating  temperature  of  the  circuit  board  was  130  °C,  greater  power  levels  were  not 
pursued. 

ET-C  was  conducted  with  three  steady  state  test  points,  1.8W,  2.7W  and  4.5W, 
which  are  the  same  with  the  previous  test.  The  circuit  board  reached  to  97.6  °C  steady 
state  temperature  at  the  4.5W  power  level.  Greater  power  levels  were  not  pursued  to 
prevent  damage.  These  data  were  used  to  determine  two  steady  state  test  points  for  the 
subsequent  tests  to  reduce  experiment  time,  because  reaching  steady  state  temperatures  in 
sample  CubeSat  experiments  required  at  least  24  hours  at  each  step. 

ET-B  and  ET-C  also  allowed  making  a  comparison  between  uncovered  and 
covered  (with  aluminum  plates)  chassis  thennal  behavior  when  using  the  same  heat 
sources  with  the  same  power  levels,  as  explained  above. 

Sample  CubeSat  Experiments  and  Simulation  (ET-D  and  CS-D) 

The  purpose  of  the  Sample  CubeSat  experiments  (ET-D)  was  to  detennine  the 
thennal  behavior  of  a  sample  CubeSat  within  the  TVAC  and  gather  enough  data  to  model 
the  CubeSat  computationally.  In  ET-D,  THCB  was  used  in  aluminum  plate  covered  6U 
chassis  (sample  CubeSat).  In  CS-D,  the  sample  CubeSat  was  modeled  in  TD  and  the 
conelated  model  was  used  to  design  a  modular  heat  dissipater.  Photographs  and  TD 
screenshots  of  the  ET-D  and  CS-D  are  shown  in  Figure  18. 


49 


(a)  The  Sample  6U  CubeSat  <b)  Thermal  Model  of  the  Sample  6U  Chassis 

includingCircuit  Board  with  Heaters  including  Circuit  Board  with  Heaters 


Figure  18:  Photographs  and  TD  Screenshots  of  the  ET-D  and  CS-D 

ET-D  includes  six  experiments  in  -15  °C  TV  AC  at  the  power  levels  1.8W  and 
4.5W  using  seven  thermocouples  in  each  test.  The  temperature  was  measured  at 
particular  points  (e.g.,  circuit  boards,  heaters,  stack  rods)  repeatedly  to  acquire 
comparable  data.  Some  uncontrolled  effects  caused  changes  at  the  temperatures  in 
successive  experiments.  For  example,  the  temperature  of  TV  AC  platen  and  TV  AC 
shroud  did  not  stay  at  the  same  level  in  every  experiment.  There  were  some  power 
fluctuations  at  the  power  source  of  the  heaters.  The  sample  CubeSat  was  disassembled 
before  each  experiment  to  place  the  thermocouples  on  different  points.  It  is  a  possibility 
that,  the  gaps  between  the  surfaces,  such  as  between  1U  stack  and  chassis  or  circuit  board 
and  1U  stack,  have  changed  in  consecutive  experiments,  and  cause  conductive  heat 
transfer  deviations  due  to  contact  resistance  changes.  This  effect  might  also  change  the 
temperatures  throughout  the  CubeSat.  As  explained  above,  the  temperature  of  particular 
points  was  measured  repeatedly  in  successive  experiments  to  acquire  dependable  data.  In 
the  case  of  acquiring  different  temperature  values  on  the  same  point,  the  most  frequently 
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occurring  temperature  was  assumed  to  be  the  reference  value  for  thermal  model  design. 
Deviated  temperatures  from  the  reference  value  on  a  particular  point  were  neglected 
taking  into  account  the  effects  expressed  above.  For  example,  the  circuit  board 
temperature  was  measured  in  five  Sample  CubeSat  4.5W  experiments.  The  results  were 
89.5  °C,  93.2  °C,  81  °C,  92.5  °C  and  93.3  °C.  It  is  evident  that  there  was  an  anomaly  in  the 
third  measurement,  so  that  outlying  value  was  ignored.  The  reference  value  was  selected 
as  93  °C  for  the  circuit  board.  All  reference  values  for  selected  points  and  their 
comparisons  to  computational  data  are  presented  in  Chapter  4. 

In  CS-D,  the  thermal  model  of  the  CubeSat  was  created  using  built-in  finite 
difference  objects.  Published  optical  properties  were  used  for  all  surfaces  except  the 
aluminum  covers  of  the  CubeSat.  The  CubeSat  covers  were  made  of  T-6061  commercial 
sheet  aluminum  with  a  published  IR  emissivity  value  of  0.1.  However,  since  some 
oxidation  occurred  on  the  covers,  the  IR  emissivity  value  of  the  covers  was  assumed  to  be 
0.15.  This  value  provided  improved  agreement  between  the  experimental  results  and 
computational  model  predictions.  For  modeling  conduction  through  chassis,  the  adjacent 
nodes  of  each  part  were  merged  to  behave  as  a  solid  body.  All  other  conductions  were 
modeled  using  the  previously  explained  the  contactor  method. 

Heaters  were  created  using  TD’s  Heater  command,  which  allows  the  user  to  apply 
a  heat  load  to  an  object  such  as  a  surface,  a  solid  or  a  single  node.  Thennoplastic 
rectangular  boxes  were  created  using  Solid  Brick,  which  are  the  same  size  as  Vishay 
Resistors.  Detailed  information  about  TVAC  computational  model  is  provided  by 
Hatzung  [Hatzung,  2014]. 
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Heat  Dissipater  Design  Process 

The  experiments  and  simulations,  ET-A  through  CS-D,  established  baseline 
conditions  for  a  hypothetical  CubeSat  without  explicit  thermal  control  hardware.  Since 
the  goal  of  this  research  is  to  create  a  modular  heat  dissipation  technique,  the  design 
should  function  over  a  large  range  of  CubeSat  missions.  To  achieve  this  goal,  a  heat  sink 
was  designed  to  cover  a  PC  104  circuit  board.  PC  104  printed  circuit  board  mounting 
scheme  is  widely  used  in  CubeSats,  because  it  leaves  a  considerable  area  for  circuitry  on 
the  board  and  has  a  simple  scheme  [Pumpkin,  2013].  Assuming  there  are  high  power 
processors  on  the  circuit  board,  the  heat  sink’s  main  purpose  is  conducting  the 
processors’  heat  to  a  radiator,  while  minimizing  the  heat  flow  into  other  regions  within 
the  CubeSat.  Since  it  was  measured  in  ET-C  experiments  that  the  temperature  of  the 
bottom  face  of  the  circuit  board  is  very  close  to  the  temperature  of  the  processor,  the  heat 
sink  was  designed  to  conduct  heat  from  the  bottom  face  of  the  heat  dissipater  too.  The 
heat  sink  connects  to  a  copper  radiator  with  soldered  copper  thermal  straps.  Since  the  IR 
emissivity  of  the  copper  is  very  low,  a  black  paint  surface  finish  was  employed  for  the 
outside  face  of  the  radiator  (3M  235  Black  Photo  Tape  was  used  to  simulate  black  paint 
in  experiments).  The  inner  faces  were  left  as  bare  copper,  because  it  was  aimed  to  radiate 
heat  out  of  the  CubeSat,  while  avoiding  radiating  heat  to  inner  surfaces  and  other 
components.  The  design  was  first  modeled  computationally.  Then  the  prototype  was  built 
and  validation  experiments  were  conducted  in  the  TVAC.  Finally,  orbit  performance 
simulations  of  the  validated  heat  dissipater  model  were  conducted  in  TD. 
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Computational  Design  of  the  Heat  Dissipater 

First,  a  copper  heat  sink  was  modeled  using  finite  difference  objects  to  cover  the 
PC  104  circuit  board.  The  first  heat  sink  design  is  shown  in  Figure  19.  Some  layers  of  the 
CubeSat  thermal  model  are  turned  off  in  each  image  to  get  a  better  view  below. 


(a)  Heat  Dissipaterwithout  Lid 


(b)  Heat  Dissipaterwith  Lid 


Topcover  (lid)  of 
the  Heat  Sink 


Radiator 

Connection 

Side 


1U  Stack 


(c)  Heat  Dissipaterwithout  Edges 


Circuit  Board 


Figure  19:  Three  Views  of  the  First  Heat  Dissipater  Design  in  the  Sample  CubeSat 

The  first  design  was  a  copper  box  sitting  in  the  1U  stack,  which  holds  the  circuit 
board.  The  heat  sink  was  designed  slightly  longer  in  the  radiator  side  to  be  in  direct 
contact  with  the  radiator.  The  connection  was  made  with  a  relatively  large  edge  of  the 
heat  dissipater,  which  is  labeled  as  radiator  connection  side  in  Figure  19.  A  thermal 
conductive  pad  was  designed  between  the  heat  sink  and  the  radiator  to  reduce  contact 
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resistance.  The  box  was  designed  as  one  piece  solid  except  the  lid  (the  top  cover  of  the 
box)  as  shown  in  Figure  20  (a). 


(a)  First  Heat  Dissipater  Design  without  Lid  (b)  First  Heat  Dissipater  Design  with  Lid 

Figure  20:  The  First  Fleat  Dissipater  Design 

Four  conductions  were  modeled  using  the  contactor  method  between  the  circuit 
board  and  the  heat  dissipater.  The  first  one  is  between  the  heaters  and  the  top  cover.  The 
second  one  is  between  the  circuit  board  and  the  bottom  face  of  the  box.  The  third  one  is 
between  the  box  and  the  radiator,  and  the  final  one  is  between  the  box  and  the  1U  stack 
rods.  For  the  first  three,  the  conductions  were  modeled  using  12  W/m-K  conductivity 
value,  which  is  a  common  value  for  thermal  conductive  pads,  which  were  modeled 
between  the  circuit  board  and  the  heat  dissipater  to  decrease  the  contact  resistance.  To 
model  conduction  between  the  box  and  the  1U  stack  rods,  plastic  was  used  as  the 
interfacing  material  because  of  the  plastic  washers  used  at  the  experiments  to  minimize 
the  heat  flow  through  chassis.  Thickness  of  the  each  side  of  the  box  was  determined 
taking  into  consideration  the  fact  that,  even  though  thennal  performance  is  getting  better, 
the  weight  of  the  box  also  increases  with  increasing  thicknesses. 
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After  the  box  dimensions  are  decided,  the  heat  dissipater  was  modeled  in 
SOLIDWORKS8  and  frequency  analysis  was  conducted  computationally  by  Mr.  Philip 
D.  Smith.  He  conducted  simulations  with  various  materials  for  mechanical  comparison. 
Although  this  baseline  geometry  was  not  particularly  problematic  in  terms  of  rigidity, 
opportunities  were  identified  for  improvement  from  the  frequency  analysis.  Since  the 
main  purpose  of  the  edges  of  the  box  is  supporting  the  faces  of  the  heat  dissipater,  instead 
of  removing  heat  from  the  heat  source,  aluminum  was  selected  for  edge  material  instead 
of  copper  to  decrease  the  weight  and  increase  the  stiffness. 

To  increase  the  manufacturability  and  provide  even  greater  rigidity,  the  material 
and  the  dimensions  of  the  box  were  redesigned.  The  dimensions  of  the  box  were  reduced 
to  the  size  of  the  circuit  board  as  shown  in  Figure  21.  The  results  of  Mr.  Smith’s 
frequency  analysis  are  shown  in  Figure  22  .  The  technical  drawing  of  the  heat  dissipater 
box  is  presented  in  Appendix  A:  Technical  Drawing  of  the  Heat  Dissipater  Box. 
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Topof  the  Heat 
Dissipater 


Aluminum  Edge  of 
the  Heat  Dissipater 


Circuit  Board 


Aluminum  Edge  of 

the  Heat  Dissipater 

Bottom  of  the  Heat 
Dissipater 

Figure  21:  The  Final  Design  of  the  Heat  Dissipater 
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Figure  22:  Frequency  Analysis  Results  of  the  Final  Heat  Dissipater  Design 
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In  the  new  design,  the  edges  were  designed  as  aluminum  due  to  its  weight  and 
stiffness  advantages  compared  to  copper,  as  explained  above.  Furthermore,  separate 
edges,  which  allow  for  various  configurations  of  the  box  with  different  thicknesses, 
increase  the  modularity  of  the  design.  Both  faces  of  the  heat  dissipater  and  the  radiator 
were  designed  as  copper.  Thicknesses  of  the  top  face,  bottom  face,  aluminum  edges  and 
the  radiator  were  set  as  1mm.  Two  conductions  were  modeled  between  the  box  and  the 
radiator.  First  conduction  was  modeled  between  the  top  face  of  the  dissipater  and  the 
radiator,  and  the  second  one  was  between  the  bottom  face  and  the  radiator.  Since  this 
design  did  not  have  a  direct  connection  between  the  heat  dissipater  box  and  the  radiator, 
thermal  straps  were  employed  between  the  radiator  and  the  both  faces  (top  and  bottom)  of 
the  box  as  shown  in  Figure  23. 


(a)  The  Heat  DissipaterThermal  Model  (b)  The  Heat  DissipaterThermal  Model  without 

Aluminum  Edges 


Bottom  Face  of  the  Heat  Dissipater 


Figure  23:  The  Heat  Dissipater  Thermal  Model 
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The  contactor  method  was  used  to  model  soldering  of  the  thennal  straps.  Four 
conductions  were  modeled  through  the  straps:  between  the  top  face  of  the  dissipater  and 
the  top  straps,  between  the  top  straps  and  the  radiator,  between  the  bottom  face  of  the 
dissipater  and  the  bottom  straps,  and  finally,  between  the  bottom  straps  and  the  radiator. 

Experimental  Validation  of  the  Design 

A  prototype  of  the  heat  dissipater  was  built  when  satisfactory  results  were 
obtained  with  the  computational  simulations.  First,  thermal  image  of  the  circuit  board 
was  taken  using  a  FLIR  thermal  imaging  camera.  The  camera  provides  a  visual  map  of 
test  article  temperatures  using  infrared  thermography,  which  can  be  explained  as 
detection  of  radiation  in  the  long-infrared  range  of  the  electromagnetic  spectrum.  The 
purpose  was  to  ensure  the  uniformity  of  heat  distribution  on  the  circuit  board.  After 
ensuring  the  uniformity,  the  experimental  validation  experiments  were  started  with  ET-E. 
The  heat  dissipater  was  placed  in  the  sample  CubeSat  with  the  THCB  heater.  The  TVAC 
was  set  to  -15  °C  and  two  power  levels  were  (1.8W  and  4.5W)  applied  to  heaters,  which 
are  the  same  conditions  used  with  ET-D.  In  the  experiments  conducted  with  THCB, 
aluminum  edges,  shown  in  Figure  21,  were  not  used  because  the  heaters  had  the  same 
thickness.  The  reason  was  to  create  a  direct  contact  between  the  heaters  and  the  top  face 
of  the  heat  dissipater.  This  adjustment  increased  the  conduction  and  simplified  the 
validation  process.  A  TGX  Ultra  Soft  Thennal  Conductive  Pad  was  used  between  the 
heaters  and  the  top  face  of  the  heat  dissipater  and  between  the  bottom  face  and  the  circuit 
board  in  all  experiments  with  THCB  as  shown  in  Figure  24. 
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(a)  Thermal  Conductive  Pad  (b)  Thermal  Conductive  Pad  Used  between  the  Resistors 
Used  between  the  Heat  and  the  ToP  Face  of  the  Heat  Dissipater 

DissipaterBottom  Face  and 
the  Circuit  Board 


Figure  24:  Thermal  Conductive  Pad  on  the  Surfaces 

ET-E  was  conducted  with  a  bare  copper  1U  size  (10  cm  x  10  cm)  radiator. 
Thermal  straps  were  not  used  in  this  experiment.  The  reason  was  to  with  prevent  a  direct 
conductive  heat  transfer  between  the  heat  dissipater  and  the  radiator,  therefore,  to  be  able 
to  observe  the  effectiveness  of  the  straps  in  the  following  experiments. 

In  ET-F,  conduction  between  the  heat  dissipater  and  the  radiator  was  created 
using  copper  thermal  straps.  The  comparison  of  ET-F  and  ET-E  provided  the 
effectiveness  of  the  thermal  straps.  The  straps  soldered  to  both  surfaces  using  a  soldering 
iron  as  shown  in  Figure  25. 
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1U  Stack 


Radiator 


CopperStraps 


Chassis 


Heat  Dissipater 
Top  Surface 


Figure  25:  The  Thermal  Straps  Soldered  to  Heat  Dissipater  and  the  Radiator 

(ET-F) 

These  thermal  straps  consisted  of  15  mm  wide  Flat  Bare  Copper  Braid,  which  is  widely 
used  in  ground  wire  applications.  Each  of  them  has  the  dimensions  of  1  mm  xl5  mm  x  70 
mm.  Four  straps  were  used  between  each  heat  dissipater  surface  and  the  radiator.  Even 
though  the  braids  are  1  mm  thick,  the  effective  thickness  of  the  braids  was  calculated  as 
0.34  mm,  because  a  12.4  cm  x  1.5  cm  braid  piece  weighs  5.1  g,  however  a  solid  copper 
with  1mm  thickness  should  weigh  16. 6g. 

In  ET-G,  3M  235  Black  Photo  Tape  was  used  to  cover  the  radiator  to  increase  its 
radiative  heat  transfer  with  TVAC.  The  purpose  was  to  determine  the  effectiveness  of  the 
surface  finish  of  the  radiator.  Photographs  of  ET-F  and  ET-G  are  shown  in  Figure  26. 
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(a)  The  Sample  CubeSat  with  1U 
Radiator  (ET-F) 


(b)  The  Sample  CubeSat  with 
Black  Photo  Tape  Covered  1U 
Radiator  (ET-G) 


Figure  26:  The  Sample  CubeSat  with  Black  Photo  Tape  Covered  1U  Radiator 


In  ET-H,  ET-G  configuration  re-conducted  with  new  soldered  straps  shown  in 
Figure  27  (b).  The  purpose  was  to  observe  the  effectiveness  of  the  strap  connection. 


(a)  Soldered  Copper  Braids  Used  in  (b)  Two  Views  of  Soldered  Copper  Braids  Used  in  ET-H 
ET-F  and  ET-G 


Figure  27:  Soldered  Copper  Braids  Used  in  ET-F,  ET-G,  and  ET-H 
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The  braids  were  soldered  using  a  soldering  iron  for  ET-F  and  ET-G.  To  increase 
the  conduction,  a  blowtorch  was  used  to  solder  the  braids  in  ET-H.  It  can  be  seen  in 
Figure  27(a)  that,  ET-F  soldering  technique  did  not  create  a  connection  between  the 
braids  and  the  dissipater  as  good  as  the  technique  in  ET-H.  ET-H  heat  dissipater 
configuration,  which  includes  the  heat  dissipater  box  without  aluminum  edges,  total  eight 
7  cm  long  copper  soldered  braids,  and  1U  size  black  tape  covered  radiator,  was  named  as 
basic  configuration  of  the  heat  dissipater  in  this  research.  Since  it  is  assumed  as  the  basic 
configuration  of  the  heat  dissipater,  ET-H  results  were  used  to  validate  computational 
model  of  the  heat  dissipater  on  TD.  Conduction  values  of  ET-H  were  used  to  model  the 
conductions  through  braids  using  the  contactor  method.  In  CS-H,  TVAC  computational 
model  was  set  to  same  conditions  with  ET-H  and  the  same  power  levels  were  applied. 

In  ET-H,  an  additional  power  level,  11W,  was  applied  to  heaters  different  from 
the  other  experimental  validation  tests.  The  purpose  was  to  determine  maximum  power 
level  to  reach  to  93  °C  on  the  circuit  board,  which  was  the  temperature  measured  without 
the  heat  dissipater  in  ET-D. 

In  ET-I,  a  3U  size  (10  cm  x  30  cm)  radiator  was  used.  The  results  were  compared 
with  ET-H  to  determine  the  effect  of  the  radiator  size.  In  ET-J,  the  3U  size  radiator  was 
covered  with  black  tape  to  increase  its  radiative  heat  transfer  within  the  TVAC.  Figure  28 
shows  the  photographs  of  the  CubeSat  used  in  the  fifth  and  the  sixth  experiments. 
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(b)  The  CubeSat  with  3U  Size  (c)  The  CubeSat  with  Black  Tape 

Radiator,  ET-I  Covered  3U  Size  Radiator,  ET-J 


Figure  28:  Photographs  of  ET-I  and  ET-J 

In  ET-K,  a  black  tape  covered  1U  size  deployable  radiator  was  used  in  addition  to 
a  1U  size  side  radiator,  as  shown  in  Figure  29,  to  determine  the  effect  of  a  deployable 
radiator.  The  ruler  in  Figure  29  is  in  inches.  Two  straps  were  soldered  to  create  the 
conduction  between  the  top  face  of  the  dissipater  and  the  deployable  radiator.  Two  braids 
were  used  between  the  top  face  and  the  side  radiator  and  four  straps  were  used  between 
the  bottom  face  and  the  side  radiator. 
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(a)  Top  View 


(b)  Bottom  View 


Braids 
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Radiator  Top  Face 

and  the  Side 
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Radiator 


Figure  29:  The  Heat  Dissipater  with  a  1U  Size  Side  Radiator  and  a  1U  Size 

Deployable  Radiator  (ET-K) 

Both  faces  of  the  deployable  radiator  and  the  outside  face  of  the  side  radiator  were 
covered  with  black  tape  to  be  able  to  make  a  comparison  with  ET-H  and  determine  the 
effect  of  the  deployable  radiator.  Two  L  shaped  copper  fixed  hinges  were  used  to  mount 
the  deployable  radiator  to  the  side  radiator.  Since  the  hinges  are  poor  conductors  in  actual 
use,  black  tape,  kapton  tape,  and  plastic  washers  were  placed  between  the  hinges  and  the 
deployable  radiator  to  reduce  conduction  between  the  deployable  radiator  and  the 
CubeSat  as  shown  in  Figure  30  (b).  Photographs  of  ET-K  are  shown  in  Figure  30. 
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(a)  Two  views  of  ET-K 


Figure  30:  Photographs  of  ET-K 


The  design  allows  the  user  to  add  or  remove  aluminum  edges,  shown  in  Figure 
21  and  Figure  31  ,  with  respect  to  circuit  board  components.  The  ruler  in  Figure  31  is  in 
inches.  The  weight  of  all  parts  (heat  dissipater  box,  radiator  and  soldered  braids)  is  360 
grams. 
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Figure  31:  The  Heat  Dissipater  with  1U  Size  Radiator 
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Since  the  thickest  components  of  the  circuit  board  were  the  heaters,  which  are  at  the  same 
height,  in  the  experiments  with  THCB,  the  aluminum  edges  were  not  mounted  as  shown 
in  Figure  25,  Figure  27,  and  Figure  28  (a).  The  top  and  bottom  faces  of  the  heat  dissipater 
were  placed  directly  on  the  circuit  board.  This  configuration  allowed  for  a  direct  contact 
and  a  good  conduction  between  the  processors  and  the  heat  dissipater’ s  top  surface. 
However,  a  direct  contact  may  not  be  available  in  many  applications.  For  example,  there 
may  be  thicker  components  than  the  main  heat  source  as  shown  in  Figure  32  (b).  In  this 
case,  aluminum  edges  of  the  heat  dissipater  must  be  used,  as  shown  in  Figure  32  (c),  and 
the  heat  must  be  conducted  from  the  high  power  processor  to  the  top  face  of  the  heat 
dissipater,  which  sits  on  the  aluminum  edges.  To  simulate  this  configuration,  two 
experiments,  ET-L  and  ET-M,  were  conducted  to  get  further  verification  for  the  design. 
SHCB,  shown  in  Figure  14  (b)  and  Figure  32(a),  was  used  in  these  experiments. 


(a)  SHCB  (The  Circuit  Board  with  a  Vishay  10W  Capacity  Resistor) 


(b)  SHCB  with  Sample 
Components 


(c)  SHCB  with  Aluminum  Edges 
ofthe  Dissipater 


Figure  32:  Photographs  of  SHCB 
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In  ET-L,  the  sample  CubeSat  was  tested  in  the  -15  °C  TV  AC  with  3.3W  power. 
The  purpose  of  the  ET-L  was  to  get  the  thermal  behavior  of  the  sample  CubeSat  with 
SHCB.  In  ET-M,  the  heat  dissipater  with  aluminum  edges  was  placed  to  cover  the  circuit 
board  as  shown  in  Figure  33,  assuming  there  are  thicker  (higher)  components  than  the 
heater  (high  power  processor)  on  the  circuit  board  and  a  direct  contact  between  the  heat 
source  and  the  heat  dissipater’s  top  surface  is  not  possible.  The  purpose  of  ET-M  was  to 
determine  the  effectiveness  of  the  heat  dissipater,  when  a  direct  contact  between  the  heat 
source  and  the  heat  dissipater  is  not  possible. 


Wt& I* 
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(a)  The  Copper  Braids  Used 
between  the  Heater  and  the  Top 
Face  of  the  Heat  Dissipater  (The 
ruleris  in  inches) 


(b)  The  Heat  Dissipater 
Covered  the  SHCB 
(Top  Face  is  Open) 


(c)  The  Braids  Are 
Placed  between  the 
Heat  Dissipater 
and  the  Single  Heater 
(Top  Face  is  Open) 


(d)  Two  views  of  the  ET-M  (Top  Face  is  Closed) 


Figure  33:  Photographs  of  ET-M 

Conduction  was  created  using  copper  braids  between  the  heater  and  the  top  face 
of  the  heat  dissipater.  The  braids  used  in  this  experiment  are  shown  in  Figure  33(a).  (The 
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ruler  in  the  picture  is  in  inches.)  The  braids  were  placed  and  pressed  between  the  single 
heater  and  the  top  face  of  the  heat  dissipater  as  shown  in  Figure  33.  In  ET-M,  TVAC  was 
set  to  -15  °C  and  3.3W  power  applied  to  heater,  which  is  the  same  condition  with  ET-L. 

In  ET-N,  SHCB  was  used  again  in  -15  °C  TVAC  with  3.3W  power.  However,  this 
time  the  basic  configuration  of  the  heat  dissipater  (without  aluminum  edges)  was  used. 
The  purpose  of  this  test  was  to  compare  results  with  ET-M  and  determine  the  heat 
transfer  effectiveness  of  the  copper  braid  pieces  used  in  ET-M. 

As  discussed  previously,  since  ET-  D  results  showed  that  the  temperature  on  the 
bottom  face  of  the  circuit  board  is  very  close  to  the  temperatures  of  the  heaters,  the  shape 
of  the  heat  dissipater  design  was  detennined  as  a  box.  The  first  reason  was  to  be  able  to 
dissipate  the  heat  using  the  bottom  face  of  the  circuit  board  in  addition  to  the  top  face. 
The  second  reason  was  to  minimize  or  eliminate  the  heat  source’s  thennal  effect  on  the 
CubeSat’s  overall  thermal  management  with  covering  the  heat  source  with  a  box  shape 
heat  sink.  In  ET-O,  only  the  top  face  of  the  heat  dissipater  was  used  with  THCB  in  -15  °C 
TVAC.  The  photographs  of  ET-0  are  shown  in  Figure  34.  The  same  two  power  levels 
with  the  other  THCB  experiments,  4.5W  and  1.8W,  were  used.  The  purpose  of  this  test 
was  to  detennine  the  effectiveness  of  the  heat  dissipater’ s  box  shape  with  a  comparison 
between  ET-0  and  ET-H. 
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Figure  34:  Photographs  of  ET-O 


Computational  Orbit  Performance  Simulations  of  the  Design 
While  it  is  able  to  simulate  the  space  environment,  a  TV  AC  cannot  reach  actual 
deep  space  vacuum  and  temperature  conditions.  The  main  purpose  of  TVAC  testing  is  to 
validate  a  thermal  design  so  that  one  may  conduct  dependable  further  computational 
experiments  such  as  orbit  performance  simulations. 

Two  orbits  were  modeled  with  TD  for  these  simulations.  The  validated  model 
created  in  CS-H  was  used  in  orbit  performance  simulations.  Both  orbits  were  modeled 
with  500  km  circular  orbits.  It  needs  to  be  emphasized  that,  the  most  important 
assumption  of  these  simulations  was  allowing  radiators  to  face  deep  space  continuously. 
The  first  orbit  had  a  0°  beta  angle  to  simulate  a  sample  cold  case  as  shown  in  Figure  35. 
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Figure  35:  Two  views  of  the  0°  Beta  Angle  Orbit,  The  Sample  Cold  Case 

The  radiator  is  on  the  3U  face  (10  cm  x  30  cm  face)  and  facing  -Z  direction  on 
this  orbit.  The  sun  is  in  +X  direction  and  the  purple  circles  indicates  the  shadow  of  the 
Earth.  In  the  second  view  in  Figure  35,  the  Earth  model  was  turned  off  to  get  a  better 
view  of  the  orbit.  Three  tests  were  conducted  with  this  orbit:  CS-P,  CS-Q,  and  CS-R.  In 
CS-P,  the  sample  CubeSat  model  was  used  without  the  heat  dissipater  and  the  radiator. 
4.5W  of  power  applied  to  heaters.  The  purpose  of  this  simulation  was  to  determine  the 
thermal  behavior  of  the  sample  CubeSat  with  THCB  in  the  sample  cold  case  orbit.  In 
CS-Q,  basic  configuration  of  the  heat  dissipater  (with  1U  size  black  tape  covered 
radiator)  was  used  to  get  the  effectiveness  of  the  heat  dissipater  in  the  cold  case  orbit. 
First,  4.5W  of  power  was  applied  to  be  able  to  make  comparison  with  CS-P  and  then,  the 
power  was  increased  to  acquire  the  same  temperature  on  the  circuit  board  with  the  CS-P 
simulation.  In  CS-R  simulation,  the  black  tape  covered  3U  size  radiator  was  used  with  the 
heat  dissipater  to  determine  the  effect  of  the  radiator  size.  Again,  4.5W  was  used  for 
comparison,  and  the  power  level  was  increased  to  acquire  the  same  temperature  on  the 
circuit  board  with  the  CS-P  simulation. 


70 


The  second  orbit  was  modeled  with  a  90°  beta  angle  for  a  sample  hot  case,  which 
is  shown  in  Figure  36.  The  radiator  is  on  the  3U  face  (10  cm  x  30  cm  face)  and  facing  -Y 
direction  in  the  Figure  36.  To  simulate  better  a  hot  case,  6U  face  (20  cm  x  30  cm  face)  of 
the  CubeSat  was  oriented  to  the  sun  direction.  The  sun  is  in  +X  direction.  The  Earth 
model  was  turned  off  in  the  second  view.  Three  simulations  were  conducted  with  this 
orbit:  CS-S,  CS-T,  and  CS-U. 


Figure  36:  Two  views  of  the  90°  Beta  Angle  Orbit,  the  Sample  Hot  Case 

In  CS-S,  the  heat  dissipater  model  and  the  radiator  were  not  used.  Thermal 
behavior  of  the  sample  CubeSat  with  THCB  in  the  sample  hot  case  orbit  was  determined 
in  CS-S  with  applying  4.5W  power  to  heaters.  In  CS-T,  4.5W  power  was  applied  to 
heaters  using  the  basic  configuration  of  the  heat  dissipater  with  the  purpose  of 
determining  the  effectiveness  of  the  heat  dissipater  on  the  sample  hot  case  orbit.  After 
4.5W  simulation,  the  power  level  was  increased  to  acquire  the  same  temperature  on  the 
circuit  board  with  CS-S  to  provide  the  improvement  in  power  capacity  with  the  same 
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maximum  temperature.  In  CS-U,  the  black  tape  covered  3U  size  radiator  was  used  to 
provide  the  effect  of  radiator  size  with  two  power  levels:  4.5W  and  22W. 

Summary 

This  chapter  provided  the  descriptions  of  the  methodologies  used  in  this  research. 
Both  experimental  and  computational  methodologies  of  the  thesis  were  introduced  and  a 
brief  experiment  procedure  was  presented  for  both  of  them.  The  test  articles  were 
presented  and  all  experiments  and  computational  simulations  conducted  in  this  research 
were  discussed  in  detail.  Finally,  the  heat  dissipater  design  process  was  introduced  and 
the  factors  that  drove  the  process  were  expressed. 
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IV.  Analysis  and  Results 


This  chapter  provides  the  data  collected  from  the  experiments  and  predictions 
computed  in  simulations.  The  summary  of  the  experiments  and  the  simulations  is 
presented  in  Table  1.  As  explained  previously,  in  this  thesis,  experimental  test  is 
abbreviated  as  ET,  computational  simulation  is  abbreviated  as  CS,  and  all  experiments 
and  simulations  were  named  with  a  letter  to  prevent  confusion.  For  example  chassis-only 
experimental  test  was  named  as  experimental  test-A  ,  ET-A,  and  chassis-only 
computational  simulation  was  named  as  computational  simulation-A,  CS-A. 

The  results  of  ET-A  and  CS-A  provided  data  for  further  validation  for  the 
previously  created  TVAC  computational  model.  ET-B  provided  data  about  THCB 
thennal  behavior  under  vacuum.  ET-C  data  were  used  to  determine  power  levels  for  the 
subsequent  experiments.  ET-D  and  CS-D  experiment  results  were  used  to  acquire  a 
validated  sample  CubeSat  model,  which  was  used  to  create  a  heat  dissipation  model. 
Lastly,  all  other  data  are  presented  in  this  chapter  provided  essential  data  for  the  objective 
of  this  thesis,  which  is  developing  a  modular  approach  to  dissipate  heat  from  heat  sources 
on  a  circuit  board  in  a  CubeSat. 

Results  of  Chassis  Only  Experiment  and  Simulation  (ET-A  and  CS-A) 

These  tests  included  three  steady  state  temperature  set  points,  -10,  30,  and  60  °C, 
which  were  selected  within  the  T VAC’s  temperature  envelop.  ET-A  data  were  used  to 
determine  the  thennal  behavior  of  the  chassis  in  TVAC.  In  CS-A,  the  chassis  were 
modeled  on  TD  and  the  results  were  used  for  validation. 

Five  thermocouples  were  used  to  collect  temperature  data  from  the  chassis.  For 

the  computational  simulations,  the  nodes  were  determined  which  are  at  the  same  location 
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as  the  thermocouples.  Selected  thermocouple  positions  are  shown  in  Figure  37. 
Temperature  color  post-processing  figures  of  the  -10  °C  CS-A  are  shown  in  Figure  38. 
Since  there  are  some  energy  losses  because  of  the  opening  to  solar  simulator  gate  (shown 
in  Figure  3),  which  is  in  the  +X  direction  in  Figure  38,  the  temperatures  are  slightly 
different  around  the  chassis.  There  are  also  some  losses  from  the  front  door,  -Y  direction, 
which  can  be  detectable  from  the  color  difference  on  the  chassis  between  the  -Y  side  and 
+Y  side.  Even  though  there  are  losses  from  the  walls  of  the  shroud,  these  losses  can  be 
assumed  to  be  compensated  since  the  walls,  except  the  front  wall,  have  active 
temperature  control  with  the  engaged  thermal  media  loop  piping.  Result  comparison  of 
the  -10  °C  ET-A  and  CS-A  is  providing  in  Table  2. 


(a)  Experimental  (b)  First  View  of  the  (c)  Second  View  of  the 

TestView  Thermocouple  Positionsfrom  Thermocouple  Positionsfrom 

the  Computational  Test  the  Computational  Test 


Figure  37:  Thermocouple  Positions  for  ET-A  and  CS-A 
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Figure  38:  Temperature  Color  Post-processing  Figures  of  the  -10  °C  CS-A 


Table  2:  T=  -10  °C  ET-A  and  CS-A  Steady  State  Temperature  Result  Comparison 


TVAC  Plate  and  Shroud  Temperature  T  =  -10  °C 

Experimental 

Results  (°C) 

Computational 

Results  (°C) 

Temperature 
Error  (°C) 

1st  Thermocouple 

-6.43 

-7.04 

0.61 

2nd  Thermocouple 

-6.45 

-7.79 

1.34 

3rd  Thermocouple 

-5.89 

-7.73 

1.84 

4th  Thermocouple 

-5.43 

-3.81 

1.62 

5th  Thermocouple 

-6.16 

-8.84 

2.68 

Table  2  shows  that,  there  is  a  good  validation  for  the  model,  with  an  error  below  2.7  °C. 

Temperature  color  post-processing  figures  of  the  30  °C  CS-A  are  shown  in  Figure 
39,  result  comparison  of  the  30  °C  ET-A  and  CS-A  is  providing  in  Table  3. 
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Figure  39:  Temperature  Color  Post-Processing  Figures  of  the  30  °C  CS-A 


Table  3:  T=  30  °C  ET-A  and  CS-A  Steady  State  Temperature  Results  Comparison 


TVAC  Plate  and  Shroud  Temperature  T  =  30  °C 

Experimental 

Results  (°C) 

Computational 
Results  (°C) 

Temperature 
Error  (°C) 

1st  Thermocouple 

30.27 

29.73 

0.54 

2nd  Thermocouple 

29.71 

29.75 

0.04 

3rd  Thermocouple 

30.46 

29.51 

0.95 

4th  Thermocouple 

30.83 

29.76 

1.07 

5th  Thermocouple 

30.03 

29.85 

0.18 

It  makes  sense  that  the  energy  losses  from  the  shroud  were  small  in  the  30  °C  test, 
which  is  close  to  room  temperature.  As  a  result,  the  error  possibility  of  the  computational 
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model  was  small  and  the  validation  for  the  model  was  very  good  with  an  error  below 
1.1  °C. 

Figure  40  shows  the  temperature  color  post-processing  figures  of  the  60  °C  ET-A 
and  CS-A.  Since  the  energy  losses  are  more  in  the  directions  of  +X  (solar  simulator  gate) 
and  -Y  (front  door),  as  expected,  the  temperatures  of  the  chassis  are  relatively  cooler  at 
those  sides.  Table  4  provides  the  result  comparison  of  the  60  °C  ET-A  and  CS-A.  Result 
comparison  of  the  60  °C  ET-A  and  CS-A  is  providing  in  Table  4. 


Figure  40:  Temperature  Color  Post-Processing  Figures  of  the  60  °C  ET-A  and  CS-A 
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Table  4:  T=  60  °C  ET-A  and  CS-A  Steady  State  Temperature  Results  Comparison 


TVAC  Plate  and  Shroud  Temperature  T  =  60  °C 

Experimental 

Results  (°C) 

Computational 
Results  (°C) 

Temperature 
Error  (°C) 

1st  Thermocouple 

58.88 

58.09 

0.79 

2nd  Thermocouple 

59.60 

58.26 

1.34 

3rd  Thermocouple 

59.33 

58.32 

1.01 

4th  Thennocouple 

58.13 

58.28 

0.14 

5th  Thennocouple 

59.71 

59.02 

0.69 

As  shown  in  Table  4,  correlation  between  the  computational  and  experimental 
60  °C  tests  was  done  with  a  temperature  error  below  1.4  °C.  Temperature  errors  for  all 
steady  state  set  points  are  provided  below  in  Table  5  for  an  easier  evaluation. 


Table  5:  Temperature  Errors  of  the  ET-A  and  CS-A  for  All  Steady  State  Test  Points 


Temperature  Errors  of  ET-A  and  CS-A  (°C) 

Steady  State  Set  Point 

T  = -15 °C 

H 

II 

w 

o 

o 

o 

T  =  60  °C 

1 st  Thermocouple 

0.61 

0.54 

0.79 

2nd  Thermocouple 

1.34 

0.05 

1.34 

3rd  Thermocouple 

1.84 

0.95 

1.01 

4th  Thermocouple 

1.62 

1.07 

0.14 

5  th  Thermocouple 

2.68 

0.18 

0.69 
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Table  5  shows  that  the  correlation  was  done  for  all  steady  state  set  points  with  an 
error  below  2.7  °C.  As  the  main  result  of  ET-A  and  CS-A,  the  greater  validation  of  the 
previously  generated  TVAC  environment  computational  model  was  provided  by 
conducting  experiments  of  a  sample  6U  CubeSat  chassis. 

Results  of  Chassis  with  Heaters  Experiment  (ET-B)  and  Sample  CubeSat  Power 
Level  Determination  Experiment  (ET-C) 

In  ET-B,  the  TVAC  was  set  to  -15  °C  and  the  experiment  was  started  with  a 
power  level  of  1.8  W.  When  steady  state  was  achieved,  the  power  was  increased  to 
2.7  W,  4.5  W,  and  6.3  W  and  data  were  collected  at  each  new  level. 

In  ET-C,  the  TVAC  was  set  to  -15  °C  and  1.8  W,  2.7  W,  and  4.5  W  power  levels 
were  applied  to  observe  the  effect  of  high  reflectivity  of  the  aluminum  covers  with 
respect  to  black  paint  coated  TVAC  shroud.  In  this  experiment,  different  from  the  chassis 
with  heaters  experiment,  powers  greater  than  4.5  W  were  not  pursued  to  prevent  damage. 

When  there  is  a  heat  source  inside  the  TVAC,  the  temperatures  of  the  shroud  and 
the  plate  change  slightly  during  long  experiments  such  as  24  hours  or  longer.  If  the 
temperature  inside  of  the  TVAC  is  relatively  large,  such  as  a  resistor  at  100  °C,  the 
change  of  the  TVAC  temperature  also  becomes  large  (±  5  °C  in  24  hours).  For  this 
reason,  since  most  of  the  experiments  in  this  research  were  longer  than  24  hours,  plate 
and  shroud  temperatures  of  the  TVAC  were  also  presented  with  thermocouple  data.  For 
example,  even  though  the  set  temperature  of  TVAC  was  -15  °C  in  ET-B,  shroud 
temperature  increased  to  -0.6  °C  when  the  steady  state  was  achieved  at  the  6.3W  level. 
Table  6  shows  the  steady  state  results  for  each  power  level. 
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Table  6:  Steady  State  Result  Comparison  of  the  ET-B  and  ET-C 


ET-B 

ET-C 

1.8W 

Power 

TV  AC 

Plate  /  Shroud  Temperature 

-14.3°C  /  -12.4°C 

-9.1°C  /  -4.3°C 

Circuit  Board  Temperature 

31.3°C 

52.3°C 

2.7W 

Power 

TV  AC 

Plate  /  Shroud  Temperature 

-10.3°C  /  -5.5°C 

-7.0°C  /  -1.9°C 

Circuit  Board  Temperature 

51.9°C 

75.0°C 

4.5W 

Power 

TV  AC 

Plate  /  Shroud  Temperature 

-7.6°C/-2.0°C 

-5.5°C  /  -0.4°C 

Circuit  Board  Temperature 

74.1°C 

102. 0°C 

6.3W 

Power 

TV  AC 

Plate  /  Shroud  Temperature 

-6.0°C  /  -0.6°C 

- 

Circuit  Board  Temperature 

92.2°C 

- 

Table  6  shows  that  temperatures  increased  faster  in  the  sample  CubeSat 
experiment  as  expected.  Resistors  in  the  aluminum-covered  chassis  heated  the  circuit 
board  to  102  °C  with  4.5W,  although  in  ET-B,  6.3W  was  needed  to  reach  to  the  92.2  °C 
temperature  level. 

Two  power  levels,  1.8W  and  4.5W,  were  selected  for  the  subsequent  experiments 
to  decrease  the  experiment  time,  taking  into  consideration  the  results  of  ET-C. 

Results  of  Sample  CubeSat  Tests  (ET-D  and  CS-D) 

Before  the  ET-D  experiments,  thermal  image  of  the  circuit  board  was  taken  to 
ensure  the  unifonnity  of  heat  distribution  on  the  circuit  board,  as  shown  in  Figure  41. 
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Figure  41:  Thermal  Image  of  the  Circuit  Board  in  Ambient  Pressure  and  Room 

Temperature 

It  should  be  noted  that,  the  circuit  board  was  perpendicular  to  the  ground.  Upper 
heaters  seem  warmer  because  of  natural  convection  effect;  however,  it  is  apparent  that, 
there  is  horizontally  uniform  heat  dissipation.  Consequently,  the  thermal  image  showed 
that  there  was  a  roughly  uniform  heat  distribution  on  the  circuit  board  and  all  resistors 
dissipate  approximately  the  same  amount  of  heat.  Therefore,  in  the  case  of  acquiring 
different  temperatures  from  two  thennocouples  on  the  circuit  board’s  top  face  in  the  same 
experiment,  assuming  uniform  heat  dissipation  on  the  board,  it  was  known  that  there  was 
an  error  with  measurements. 

Six  ET-D  experiments  were  conducted  with  1.8W  and  4.5W  power  levels.  In 

CS-D,  the  same  experiments  were  conducted  in  TD  as  well.  Steady  state  result 

comparison  for  1.8W  ET-D  and  CS-D  is  shown  in  Figure  42.  Steady  state  result 
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comparison  for  4.5W  is  shown  in  Figure  43.  Some  layers  of  the  CubeSat  thennal  model 
were  turned  off  to  get  a  better  view  of  the  salient  features  in  Figure  42  and  Figure  43. 
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Twctrature  1 Q3  ^  T'ne  =0  sec. 


Chassis  Middle: 
ER:  20. CPC 
CR:  22.0 'C 


1U  Stack  Top  (+X): 
ER:  34.5°C 
CR:  31.0°C 


(b)  Inside  View-1 


Circuit  Board 
Top  Face: 
ER:  93.0  'C 
CR:  90.0°C 


Spacer: 

ER:  48.0°C 
CR:  44.0°C 


(c)  Inside  View-2 


Circuit  Board 
Back  Face  : 
ER: 92 CC 
CR:  89.5:C 


Heater: 

ER:  94.0°C 
CR:  92.0:’C 


Figure  43:  Steady  State  Test  Result  Comparisons  for  the  4.5W  ET-D  and  CS-D 

84 


Table  7  shows  the  result  comparison  for  1.8W  and  Table  8  for  4.5W  ET-D  and 
CS-D.  Each  measurement  point,  e.g.  front  cover,  back  cover,  or  chassis  middle,  is 
represented  with  a  number  in  the  tables. 


Table  7:  1.8W  ET-D  and  CS-D  Steady  State  Result  Comparison 


Power  Level  P=1.8W 

TVAC  Plate  Temperature  T  =  -14  °C,  TVAC  Shroud  Temperature  T=  -13°C 

Measurement  Point 

Experimental  Results 
(°C) 

Computational 
Results  (°C) 

Error  (°C) 

Measurement  Point 

Experimental  Results 
(°C) 

Computational 
Results  (°C) 

U 

O 

u 

© 

u 

u 

w 

1 

-0.2 

-2.1 

wm 

7 

3.5 

0.3 

3.2 

2 

-1.6 

-1.0 

0.6 

8 

9.0 

7.0 

2.0 

3 

2.2 

0.1 

wm 

9 

8.0 

5.5 

2.5 

4 

0.2 

-1.5 

1.7 

10 

37.0 

39.0 

2.0 

5 

5.0 

2.7 

2.3 

11 

15.4 

14.0 

H 

6 

38.0 

40.0 

2.0 

12 

36.0 

37.5 
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Table  8:  4.5W  ET-D  and  CS-D  Steady  State  Result  Comparison 


Power  Level  P=4.5W 

TVAC  Plate  Temperature  T  =  -15  °C,  TVAC  Shroud  Temperature  r 

r=  -i4°c 

Measurement  Point 

Experimental  Results 
(°C) 

Computational 

Results  (°C) 

Error  (°C) 

Measurement  Point 

Experimental  Results 
(°C) 

Computational 

Results  (°C) 

Error  (°C) 

1 

13.0 

10.3 

2.7 

7 

90.0 

93.0 

3.0 

2 

11.2 

13.0 

1.8 

8 

89.5 

92.0 

2.5 

3 

16.0 

13.8 

2.2 

9 

48.0 

44.0 

4.0 

4 

14.0 

11.1 

2.9 

10 

92.0 

94.0 

2.0 

5 

22.0 

20.0 

2.0 

11 

22.0 

19.6 

2.4 

6 

34.5 

31.0 

3.5 

Table  7  and  Table  8  show  that  maximum  steady  state  temperature  error  between 
ET-D  and  CS-D  was  Terror  <  4  °C,  which  is  within  the  acceptable  error  value  determined 
for  this  research.  Consequently,  physical  shape  and  thermo-optic  properties  of  the 
computational  model  were  close  enough  to  physical  tests  and  the  computational  model 
was  dependable  for  a  heat  dissipater  design. 

Experimental  and  Computational  Results  of  the  Heat  Dissipater  Design 

Steady  state  temperature  results  of  the  heat  dissipater  validation  experiments  with 
the  THCB  are  shown  in  Table  9  for  4.5W  and  1.8W  power  levels. 
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Table  9:  Steady  State  Temperature  Results  of  the  Heat  Dissipater  Validation 
Experiments  with  THCB 


Power  Level  P=4.5W 


Thermocouple 

ET-E 

(°C) 

ET-F 

(°C) 

ET-G 

(°C) 

ET-H 

(°C) 

ET-I 

(°C) 

ET-J 

(°C) 

ET-K 

(°C) 

ET-0 

(°C) 

ET-D 

(°C) 

Circuit  Board 
Top  Face  (Two 
Thermocouples) 

103.1/ 

105.4 

56.1/ 

57.4 

46.5/ 

47.8 

34.4/ 

34.9 

43.9/ 

44.4 

26.4/ 

27.3 

30.0/ 

31.5 

69.5/ 

67.5 

93.0 

Heat  Dissipater 
Top  Face 

(Failed) 

46.4 

36.9 

33.7 

39.2 

22.0 

23.3 

42.5 

- 

Heat  Dissipater 
Bottom  Face 

90.8 

49.3 

40.0 

31.8 

36.7 

19.3 

28.6 

- 

- 

Radiator 
(Outside  Face) 

22.1 

23.3 

16.4 

18.1 

26.6 

9.3 

10.1/3.8 

(Deployed) 

14.8 

14.0 

(Aluminum 

cover) 

TVAC  Plate 

-15.1 

-15.3 

-15.3 

-15.3 

-14.1 

-15.3 

-14.8 

-13.5 

-15.0 

TV  AC  Shroud 

-14.7 

-15.0 

-13.0 

-15.0 

-15.1 

-14.9 

-13.2 

-13.0 

-14.0 

Power  Level  P=1.8W 

Thermocouple 

ET-E 

(°C) 

ET-F 

(°Q 

ET-G 

(°C) 

ET-H 

(°C) 

ET-I 

(°C) 

ET-J 

(°C) 

ET-K 

(°C) 

ET-0 

(°C) 

ET-D 

(°C) 

Circuit  Board 
Top  Face  (Two 
Thermocouples) 

46.5/ 

47.5 

18.9/ 

19.5 

14.0/ 

14.6 

13.7/ 

13.9 

12.8/ 

13.2 

3.7/ 

4.1 

9.1/ 

9.8 

27.4/ 

28.6 

39.0 

Heat  Dissipater 
Top  Face 

(Failed) 

15.0 

10.2 

13.7 

11.2 

2.3 

6.5 

16.2 

- 

Heat  Dissipater 
Bottom  Face 

41.6 

16.3 

11.6 

12.6 

10.3 

1.3 

8.5 

- 

- 

Radiator 
(Outside  Face) 

4.9 

5.1 

1.9 

7.4 

6.5 

-2.2 

1.4/  -2.2 
(Deployed) 

3.7 

-1.5 

(Aluminum 

cover) 

TVAC  Plate 

-13.4 

-14.9 

-14.6 

-9.3 

-13.9 

-15.1 

-13.0 

-12.0 

-14.0 

TVAC  Shroud 

-12.9 

-14.0 

-13.8 

-5.6 

-12.6 

-14.3 

-11.8 

-11.5 

-13.0 
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Since  both  power  level  experiments  resulted  in  similar  consequences,  only  4.5W 
results  are  discussed  in  this  chapter.  ET-E  results  showed  that,  as  expected,  the  circuit 
board  temperature  was  higher  with  respect  to  ET-D,  because  the  black  heaters’  view  to 
aluminum  of  the  CubeSat  was  blocked  by  bare  copper,  top  face  of  the  heat  dissipater, 
which  has  lower  infrared  emissivity  and  absorptivity.  In  addition  to  that,  a  1U  size  bare 
copper  radiator  was  placed  instead  of  a  same  size  aluminum  plate,  as  shown  in  Figure 
26  (a),  which  has  lower  IR  emissivity  than  the  aluminum.  Since  the  conduction  between 
the  heat  dissipater  and  the  radiator  was  not  created  using  braids  in  this  experiment,  as 
explained  in  the  previous  chapter,  the  circuit  board  temperature  was  increased  to  105  °C 
from  the  93  °C  (ET-D)  at  the  4.5W  power  level. 

In  ET-F,  copper  braids  were  soldered  between  the  heat  dissipater  and  the  1U  size 
radiator  as  shown  in  Figure  25.  The  conduction  through  copper  braids  reduced  the  circuit 
board  temperature  from  105.4  °C  to  57.4  °C.  Since  copper  has  a  very  low  IR  emissivity 
(£Copper=0.04),  in  ET-G,  the  radiator  was  covered  with  black  tape  which  has  a  much  higher 
IR  emissivity  (ebiack  tape=0.9).  This  configuration  reduced  the  temperature  an  additional 
10  °C,  from  57.4  °C  to  47.8  °C.  It  must  be  noted  that,  in  ET-F  and  ET-G,  the  top  face  of 
the  heat  dissipater  is  cooler  than  the  bottom  face,  which  indicates  that  there  is  not  good 
conduction  between  the  bottom  face  of  the  heat  dissipater  and  the  radiator.  Since  the 
heaters  were  hotter  than  the  bottom  face  of  the  circuit  board  in  ET-D,  it  should  be 
expected  that  the  top  face  of  the  heat  dissipater  must  be  warmer  than  the  bottom  face, 
because  the  heaters  were  in  direct  contact  with  the  top  face  using  a  thermal  conductive 
pad.  Thus,  it  can  be  deduced  that,  there  was  a  poor  conduction  between  the  bottom  face 
of  the  heat  dissipater  and  the  radiator. 
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In  ET-H,  braids  are  soldered  at  a  very  high  temperature  using  a  blowtorch  instead 
of  soldering  iron  to  create  better  conduction.  Black  tape  covered  1U  size  radiator  was 
used,  which  is  the  same  configuration  with  ET-G.  Results  of  ET-H  showed  that  this 
soldering  technique  was  better,  as  expected.  The  temperature  difference  between  the  top 
and  bottom  faces  of  the  heat  dissipater  was  reduced  to  1.9  °C  from  3.1  °C,  which  shows 
that  there  is  very  good  conduction  between  the  radiator  and  both  faces  of  the  heat 
dissipater  through  copper  braids.  Additionally,  the  temperature  difference  between  the 
heat  dissipater  and  the  radiator  was  reduced  to  13.7  °C  from  ET-G  result  of  20.5  °C, 
which  is  another  indicator  of  better  conduction  through  braids.  Furthennore,  in  ET-G,  the 
top  face  of  the  heat  dissipater  was  hotter  than  the  bottom  face,  which  shows  that  the 
conduction  from  both  faces  of  the  heat  dissipater  to  radiator  was  relatively  equal. 
Consequently,  ET-H  results  showed  that  the  connection  technique  of  the  thennal  straps  is 
very  effective  to  determine  the  achievement  of  the  heat  dissipation.  ET-H  heat  dissipater 
configuration  was  named  as  basic  configuration  as  explained  in  the  previous  chapter. 

As  the  last  step  in  ET-H,  the  power  level  was  increased  to  acquire  a  circuit  board 
top  face  temperature  of  93  °C  (ET-D  result).  The  purpose  was  to  detennine  how  much 
power  increase  could  be  achieved  by  using  the  heat  dissipater,  assuming  the  ET-D  result, 
93  °C,  is  the  maximum  operating  temperature. 
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Table  10:  The  Comparison  of  ET-D  and  ET-H  for  Different  Power  Levels 


ET-D 

ET-H 

ET-H 

Power  (W) 

4.5 

4.5 

11 

Circuit  Board  Top  Face  Temperature 
(°C) 

93.0 

34.9 

89.0 

TVAC  Plate  Temperature  (°C) 

-15.0 

-15.3 

-13.2 

TVAC  Shroud  Temperature  (°C) 

-14.0 

-15.0 

-11.2 

Table  10  shows  that,  THCB  power  level  could  be  increased  to  1 1W  using  the  heat 
dissipater  assuming  the  maximum  temperature  allowed  is  93  °C.  It  should  be  noted  that, 
in  11W  ET-H  experiment  TVAC  was  wanner  then  ET-D  experiment,  which  indicates 
that,  actual  power  increase  capability  of  the  heat  dissipater  is  a  more  than  1 1W  for  basic 
configuration.  Consequently,  the  heat  dissipater  basic  configuration  is  capable  to  allow  an 
increase  of  power  more  than  240%  on  the  circuit  board,  without  exceeding  the 
temperature  level  without  a  heat  dissipater. 

In  ET-I,  3U  size  bare  copper  radiator  made  cooler  all  measured  temperatures  on 
the  CubeSat  with  respect  to  ET-F  as  shown  in  Table  9.  In  ET-J,  the  3U  size  radiator  was 
covered  with  black  tape  and  the  circuit  board  temperature  was  measured  as  27.3  °C, 
which  was  93.0  °C  in  ET-D  without  the  heat  dissipater. 

In  ET-K,  a  1U  size  deployable  radiator  was  added  to  basic  configuration,  and  the 
circuit  board  temperature  was  measured  as  31.5  °C,  which  is  a  better  result  with  respect 
to  ET-H  result  of  34.9  °C.  The  top  face  of  the  heat  dissipater  was  colder  than  the  bottom 
face,  which  shows  the  braids  between  the  deployable  radiator  and  the  top  face  of  the  heat 
dissipater  achieved  their  design  goal.  However,  it  is  obvious  that  the  deployable  radiator 
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did  not  improve  significantly  the  heat  dissipation  capability  of  the  design.  Even  though 
the  side  radiator  was  10.1  °C,  the  deployable  radiator  was  3.8  °C,  which  shows  that  there 
is  not  enough  conduction  between  the  top  surface  of  the  heat  dissipater  and  the 
deployable  radiator.  The  conduction  between  the  deployable  radiator  and  the  side 
radiator  was  minimized  using  black  tape,  kapton  tape,  and  plastic  washers  as  shown  in 
Figure  30  (b)  to  simulate  a  real  hinge  as  discussed  in  the  previous  chapter.  It  can  be 
deduced  that  the  effectiveness  of  the  design  with  deployable  radiator  can  be  improved  by 
increasing  the  conduction  between  the  deployable  radiator  and  the  heat  dissipater.  It 
should  also  be  noted  that  in  ET-K,  the  shroud  temperature  is  1.8  °C  warmer  than  ET-H, 
which  is  another  reason  for  not  having  a  noteworthy  result  from  ET-K. 


Table  11:  Steady  State  Temperature  Results  of  ET-L,  ET-M  and  ET-N 


ET-L 

ET-M 

ET-N 

Heat  Dissipater 

- 

V 

V  (w/o  Edges) 

Power 

3.3W 

3.3W 

3.3W 

Circuit  Board 
Temperature 

100.7  °C 

37.7  °C 

25.2  °C 

TVAC  Plate 
Temperature 

-13.7  °C 

-14.7  °C 

-14.9  °C 

TVAC  Shroud 
Temperature 

-12.9  °C 

-13.4  °C 

-13.4  °C 

Steady  state  temperature  results  of  ET-L  and  ET-M  are  presented  in  Table  11.  In 
ET-L  and  ET-M,  the  SHCB  was  used  to  determine  the  effectiveness  of  the  design  when  a 
direct  connection  between  the  heat  source  and  the  top  face  of  the  heat  dissipater  is  not 
possible,  as  discussed  in  the  previous  chapter.  First,  in  ET-L,  the  CubeSat  with  SHCB 
was  tested  in  -15  °C  TVAC.  The  experiment  was  terminated  at  the  3.3W  power  level, 


when  the  circuit  board  reached  to  100.7  °C  steady  state  temperature  to  prevent  damage  to 
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circuitry.  In  ET-M,  aluminum  edges  were  added  to  basic  configuration  of  the  heat 
dissipater  and  braid  pieces  were  squeezed  between  the  heater  and  the  top  face  of  the  heat 
dissipater  as  shown  in  Figure  33.  The  temperature  of  the  circuit  board  was  reduced  to 
37.7  °C  with  an  indirect  contact.  This  result  shows  that,  the  heat  dissipater  design  is  very 
effective  to  reduce  the  heat  by  an  indirect  contact  between  the  heat  source  and  the  heat 
dissipater. 

In  ET-N,  basic  configuration  of  the  heat  dissipater  was  used  and  top  face  of  the 
heat  dissipater  was  placed  directly  on  top  of  the  single  heater  using  a  thermal  conductive 
pad.  The  purpose  was  to  detennine  the  heat  transfer  effectiveness  of  the  basic 
configuration  in  case  of  a  concentrated  thermal  load.  Table  1  lshows  that  circuit  board’s 
temperature  was  measured  as  25.2  °C  in  ET-N  and  it  can  be  deduced  that  the  basic 
configuration  of  the  heat  dissipater  is  very  effective  on  concentrated  thermal  loads. 

In  ET-O,  only  the  top  face  of  the  heat  dissipater  was  used  with  THCB.  The 
bottom  face  of  the  heat  dissipater  was  removed  from  the  basic  configuration.  The  purpose 
was  to  detennine  the  effectiveness  of  the  heat  dissipater’ s  box  shape  with  a  comparison 
with  ET-H.  As  discussed  previously,  the  basic  configuration  of  the  heat  dissipater 
reduced  the  circuit  board  temperature  to  34.9  °C  from  93.0  °C.  However,  temperature 
decreased  to  69.5  °C  in  ET-O.  Since  the  bottom  face  of  the  heat  dissipater  caused  an 
additional  34.6  °C  decrease  in  ET-H,  it  can  be  deduced  that  the  box  shape  is  very 
effective.  The  weight  of  one  face  of  the  heat  dissipater  with  copper  braids  is  90  grams  (a 
10  cm  x  10  cm  x  0.1  cm  copper  plate,  bottom  face,  and  four  1.5  cm  x  7  cm  x  0.1  cm 
copper  braids).  Thus,  it  can  be  concluded  that  the  increase  in  weight  was  worth  the 
increase  in  capability  of  dissipating  heat. 
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Since  the  ET-H  design  was  assumed  as  the  basic  configuration  with  a  black  tape 
covered  1U  size  radiator,  the  thennal  desktop  model  of  the  heat  dissipater  (CS-H)  was 
adjusted  to  fit  the  ET-H  data.  The  thermal  braids  were  simulated  as  a  one  piece  solid  as 
shown  in  Figure  23  with  the  same  length  and  the  width  of  the  real  braids.  Since  four 
braids  with  1.5  cm  width  used  at  the  both  faces,  the  braid  model  was  created  as  one  piece 
with  6  cm  width.  The  thickness  was  reduced  to  simulate  the  effective  thickness  of  the 
braids  for  conduction.  Conductions  were  modeled  using  the  contactor  method  with  ray 
trace  algorithm  between  the  braids-the  radiator  and  the  braids-the  heat  dissipater.  A  ray 
trace  algorithm  shoots  rays  perpendicular  from  the  “from  surface”  and  deems  the  closest 
node  to  be  whatever  node  ray  first  encounters  as  explained  in  the  previous  chapter.  This 
algorithm  is  very  fast  and  accurate  between  parallel  surfaces.  In  view  of  the  fact  that  the 
contactor  calculates  all  surface  area  between  the  parallel  surfaces  for  conduction,  which 
is  not  realistic  for  this  case,  the  conductivity  values  were  determined  less  with  respect  to  a 
contactor  with  solder  intervening  material.  In  ET-H,  the  soldered  surface  areas,  which 
should  be  calculated  for  conduction,  were  smaller  as  shown  in  Figure  27  (b).  Figure  44 
shows  the  temperature  color  post-processing  figures  of  CS-H  and  Table  12  presents  the 
steady  state  temperature  result  comparison  of  ET-H  and  CS-H. 
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a)  Two  Views  of  4.5W  CS-H 


(a)  Two  Views  of  1.8W  CS-H 

Figure  44:  The  Temperature  Color  Post-Processing  Figures  of  CS-H 


Table  12:  Steady  State  Temperature  Result  Comparison  of  ET-H  and  CS-H 


ET-H  at  4.5W 
(°C) 

CS-H  at  4.5W 
(°C) 

ET-H  at  1.8W 
(°C) 

CS-H  at  1.8W 
(°C) 

Circuit  Board 
Top  Face 

34.4/ 

34.9 

34.5 

13.7/ 

13.9 

14.4 

Heat 

Dissipater  Top 
Face 

33.7 

33.6 

13.7 

13.8 

Heat 

Dissipater 
Bottom  Face 

31.8 

30.5 

12.6 

13.2 

Radiator 
(Outside  Face) 

18.1 

21.1 

7.4 

9.6 

TV  AC  Plate 

-15.3 

-15.3 

-9.3 

-9.3 

TV  AC  Shroud 

-15.0 

-15.0 

-5.6 

-5.6 
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Table  12  shows  that  the  heat  dissipater  model  was  validated  with  an  error  below 
3.0  °C  at  all  power  levels.  This  indicates  that,  the  model  is  dependable  to  conduct  further 
computational  experiments  such  as  orbit  perfonnance  simulations. 

Orbit  Performance  Simulation  Results  of  the  Heat  Dissipater 

Two  orbits  were  created  and  six  orbital  simulations  were  conducted  on  TD  as 
explained  in  the  previous  chapter.  Figure  45  shows  the  temperature  color  post-processing 
figures  of  CS-P,  CS-Q,  and  CS-R  and  the  circuit  board  steady  state  temperature 
predictions  of  the  simulations  are  presented  in  Table  13. 
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Temperature  (*C)  Temperature  (°C) 


(a)CS-P,  4.5W 


(b)CS-Q,  4.5W 


(c)  CS-Q,  13.5W 


(e)  CS-R,  19W 


(d)CS-R,  4.5W 


Figure  45:  Temperature  Color  Post-Processing  Figures  of  CS-P,  CS-Q,  and  CS-R 
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Table  13:  Steady  State  Temperature  Predictions  of  CS-P,  CS-Q,  and  CS-R 


CS-P 

CS-Q 

CS-R 

Beta  Angle 

0° 

0° 

0° 

Heat  Dissipater 

- 

V 

V 

Radiator 

Size 

- 

1U  (Black  Tape  Covered) 

3U  (Black  Tape  Covered) 

Power 

4.5W 

4.5W 

13.5W 

4.5W 

19W 

Circuit  Board 
Temperature 

96.5  °C 

22.4  °C 

96.4  °C 

4.5  °C 

97.5  °C 

The  temperature  of  the  circuit  board  was  predicted  as  96.5  °C  with  4.5W  power 
on  THCB,  in  the  sample  cold  case  orbit.  In  CS-Q,  validated  heat  dissipater  model  was 
used  and  the  circuit  board  temperature  was  decreased  to  22.4  °C  with  1U  size  black  tape 
covered  (black  painted  simulation)  radiator.  Assuming  the  96.5  °C  is  the  maximum 
operating  temperature  for  this  circuit  board,  the  power  can  be  increased  to  13.5W  using 
the  heat  dissipater,  which  means  300%  increase  in  power  capacity  for  this  configuration. 
CS-P  results  shows  that,  using  a  3U  size  black  radiator,  the  circuit  board  temperature  was 
reduced  to  4.5  °C  and  the  power  can  be  increased  to  19W  without  exceeding  96.5  °C. 
CS-P  power  result,  19W,  is  420%  of  the  original  power.  Figure  46  shows  the  temperature 
color  post-processing  figures  of  CS-S,  CS-T,  and  CS-U  and  the  circuit  board  steady  state 
temperature  predictions  of  the  simulations  are  presented  in  Table  14. 
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Figure  46:  Temperature  Color  Post-Processing  Figures  of  CS-S,  CS-T,  and  CS-U 
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Table  14:  Steady  State  Temperature  Predictions  of  CS-S,  CS-T,  and  CS-U 


CS-S 

CS-T 

CS-U 

Beta  Angle 

o 

o 

90° 

90° 

Heat  Dissipater 

- 

V 

V 

Radiator 

Size 

- 

1U  (Black  Tape  Covered) 

3U  (Black  Tape  Covered) 

Power 

4.5W 

4.5W 

15. 5W 

4.5W 

22W 

Circuit  Board 
Temperature 

1 16.8  °C 

34.6  °C 

115.3  °C 

12.2°C 

1 15.7°C 

Three  simulations  were  created  with  a  sample  hot  case  orbit.  Circuit  board 
temperature  reached  to  116.8  °C  in  this  orbit  without  the  heat  dissipater.  In  CS-T,  the 
circuit  board  temperature  was  reduced  to  34.6  °C  with  the  same  power  level  of  4.5W. 
Assuming  the  maximum  operating  temperature  is  1 16.8  °C,  the  power  on  the  processors 
can  be  increased  by  340%  to  15.5W.  In  CS-U,  the  heat  dissipater  with  3U  size  radiator 
resulted  in  12.2  °C  temperature  on  the  circuit  board.  Lastly,  CS-U  results  showed  that, 
assuming  116.8°C  is  the  maximum  operating  temperature  in  this  case,  power  on  the 
processors  can  be  increased  almost  490%  to  22W  level  without  exceeding  the  assumed 
maximum  operating  temperature. 

Summary 

All  data  collected  from  the  simulations  and  experiments  were  presented  and  the 
results  were  analyzed  in  this  chapter.  ET-A  and  CS-A  results  provided  further  validation 
to  previously  generated  TVAC  computational  model.  ET-D  and  CS-D  results  validated 
the  sample  CubeSat  computational  model.  ET-E,  F,  G,  H,  I,  J,  K,  L,  M,  and  N  results 
showed  that  the  heat  dissipater  design  is  very  effective  in  dissipating  heat  from  various 
heat  sources  and  it  is  possible  to  use  the  heat  dissipater  effectively  with  different 
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configurations,  such  as  with  various  radiator  sizes  and  with/without  the  aluminum  edges. 
ET-0  data  showed  that,  the  box  shape  heat  sink  design  is  much  more  effective  than  a 
design  with  only  one  plate  on  the  heat  source.  ET-H  and  CS-H  results  validated  the  heat 
dissipater  computational  model.  Finally,  CS-P,  Q,  R,  S,  T,  and  U  simulations  proved  that 
this  heat  dissipater  gives  valuable  results  on  orbit,  which  is  also  a  proof  of  the 
achievement  of  the  objective  of  the  thesis. 
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V.  Conclusions  and  Recommendations 


This  chapter  presents  the  conclusions  of  the  research.  The  objective  is  reviewed, 
demonstrating  that  it  has  been  met.  Advantages  and  weak  points  of  the  heat  dissipater 
design  are  discussed.  Finally,  the  recommendations  for  future  work  are  presented. 

Conclusions  of  Research 

The  objective  of  this  thesis  was  to  develop  a  modular  approach  to  dissipate  heat 
from  high  power  heat  sources  on  a  circuit  board  by  using  thermal  control  hardware  in  a 
sample  CubeSat. 

Since  the  previously  generated  TVAC  environment  computational  model  was 
validated  only  by  using  a  10  in  x  10  in  aluminum  plate  [Hatzung,  2014],  this  research 
provided  data  by  conducting  experiments  and  computational  simulations  with  a  more 
realistic  test  article,  a  sample  6U  CubeSat  chassis,  (ET-A  and  CS-A)  within  TVAC.  The 
steady  state  test  points  covered  a  large  range  of  temperatures  (-10  °C  to  60  °C)  within 
TVAC’s  temperature  envelope.  Experiment  and  simulation  results  were  validated  with  a 
maximum  steady  state  temperature  error  below  2.7  °C  for  all  steady  state  test  points. 
Thus,  further  validation  for  the  TVAC  model  was  provided. 

Later,  a  circuit  board  with  ten  resistors  (THCB)  was  placed  in  chassis  and  an 
experiment  conducted  in  -15  °C  TVAC  (ET-B).  Then,  the  same  experiment  was 
re-conducted  with  a  sample  CubeSat  (ET-C).  The  chassis  was  covered  with  aluminum 
plates  to  simulate  a  fully  enclosed  sample  CubeSat  with  heat  sources.  These  two 
experiments  showed  that,  expectedly,  the  resistors  became  much  hotter  when  their  open 
view  to  highly  absorptive  black  painted  TVAC  shroud  was  obstructed  by  highly 


reflective  aluminum  plates.  For  example,  even  though  the  temperature  of  the  top  face  of 
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the  circuit  board  was  measured  as  92.2  °C  at  6.3W  power  level  in  ET-B,  it  reached  to 
102.0  °C  at  4.5W  power  in  ET-C. 

IR  images  of  the  THCB  in  ambient  pressure  showed  that  there  was  a  roughly 
uniform  heat  distribution  on  the  circuit  board  face,  so  two  thermocouples  were  used  in 
many  experiments  to  acquire  dependable  temperature  measurements  of  the  circuit  board 
face.  The  circuit  board  top  face  temperature  was  assumed  as  reference  value  in  many 
experiments  and  simulations.  The  first  reason  was  the  difficulties  of  taking  temperature 
measurements  from  the  resistors,  because  of  their  sizes  and  shapes.  The  second  reason, 
experiments  showed  that  the  resistor  and  the  circuit  board  temperatures  were  very  close 
to  each  other  in  steady  state  results.  The  sample  CubeSat  with  THCB  was  modeled  and 
validated  by  data  of  ET-D  and  CS-D,  with  a  maximum  steady  state  temperature  error 
below  4  °C  for  two  different  power  levels.  These  results  showed  that  the  sample  CubeSat 
model  was  dependable  to  model  a  heat  dissipater  design. 

Since  the  experiment  results  showed  that  the  circuit  board  bottom  face’s 
temperature  is  very  close  to  the  top  face,  the  heat  dissipater  was  designed  as  a  box  shape 
to  remove  the  heat  using  both  faces.  Additionally,  this  heat  dissipater  was  designed  to 
minimize  or  eliminate  the  high  power  heat  sources  thermal  effect  on  the  satellite’s  overall 
thennal  management  by  preventing  heat  flow  from  the  bottom  face  of  the  circuit  board  to 
other  components  of  the  CubeSat.  Furthermore,  ET-0  results  proved  that,  removing  heat 
from  a  circuit  board  using  both  faces  is  much  more  effective  than  doing  it  with  single 
face  as  expected. 

The  heat  dissipater’s  top  and  bottom  faces  were  made  of  copper,  because  of  the 
advantage  of  thermal  conductivity.  They  were  designed  as  10  cm  x  10  cm  x  1  cm  plates. 
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Even  though  computational  simulation  results  for  the  thicker  plates  were  relatively  more 
effective  at  dissipating  heat,  weight  and  mechanical  strength  were  other  factors  that 
should  be  considered.  When  the  copper  plates  were  designed  thicker,  the  weight  became 
much  higher  because  of  the  high  density  of  copper.  Since  the  main  purpose  of  the  edges 
is  supporting  the  faces  of  the  heat  dissipater,  instead  of  removing  heat  from  the  heat 
source,  they  were  designed  as  aluminum  to  take  advantage  of  the  aluminum’s  low  density 
and  superior  mechanical  strength  with  respect  to  copper. 

The  heat  dissipater  box  was  connected  to  a  radiator  using  copper  braids,  which  are 
commonly  used  for  ground  wire  applications.  In  this  research,  four  braids  with  1.5  cm 
width  were  used  between  each  face  of  the  heat  dissipater  box  and  the  radiator.  Braid 
amount  can  be  adjusted  for  various  cases  with  different  power  levels.  Since  many 
experiments  were  conducted  with  the  heat  dissipater,  braids  were  designed  relatively  long 
(7  cm  each)  to  simplify  assembly.  Since  the  soldered  parts  were  no  longer  than  2  cm  at 
each  side  (soldered  parts  between  braids-radiator  and  braids-heat  dissipater  faces),  this 
length  was  obviously  a  disadvantage  for  the  heat  transfer. 

Even  though  the  prototype  was  built  using  commercial  copper  plates  and  copper 
braids  using  for  ground  wire  applications,  and  the  computational  models  were  designed 
based  on  these  material  properties,  the  experiment  results  showed  that,  the  heat  dissipater 
is  very  effective  to  dissipate  heat  from  a  circuit  board  containing  heat  sources  in  a 
CubeSat.  For  example,  the  design  reduced  the  temperature  of  a  circuit  board  with  ten 
resistors  from  93  °C  to  34.4  °C  at  the  4.5W  power  level,  with  1U  size  (10  cm  x  10  cm) 
radiator,  in  -15  °C  TV  AC.  In  addition  to  this,  it  reduced  the  temperature  of  a  circuit  board 
with  a  single  resistor  simulating  a  concentrated  heat  load,  with  3.3W  applied  of  power, 
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from  100.7  °C  to  25.2  °C,  with  1U  size  radiator,  in  -15  °C  TVAC.  Undoubtedly,  using 
higher  thermally  conductive  materials,  such  as  OFHC  copper  or  graphite  composites,  can 
improve  the  results. 

The  connection  technique  of  the  thermal  straps  is  also  very  important,  as  the 
comparison  of  ET-G  and  ET-H  results  showed.  In  the  first  place,  thermal  straps  were 
soldered  using  a  soldering  iron.  The  circuit  board  temperature  was  reduced  from  93  °C  to 

47.8  °C  using  the  heat  dissipater  in  a  test  configuration.  However,  when  the  braids  were 
soldered  with  high  temperature,  using  a  blowtorch,  the  temperature  was  measured  as 

34.9  °C  with  the  same  test  configuration,  which  is  a  proof  of  the  importance  of  the 
connection  quality  of  the  thermal  straps. 

The  design  is  also  very  effective  when  a  direct  contact  between  a  face  of  the  heat 
dissipater  and  the  heat  source  is  not  possible.  Using  some  thermally  conductive  materials, 
such  as  squeezed  copper  braid  pieces,  gave  very  good  results.  For  example,  first,  a  test 
conducted  with  a  single  resistor  gave  100.7  °C  on  the  circuit  board  without  the  heat 
dissipater.  Then  the  heat  dissipater  reduced  the  temperature  to  37.7  °C  in  the  same  test 
configuration,  with  an  indirect  contact  between  the  resistor  and  the  heat  dissipater  box. 
Since  the  heat  dissipater  gave  the  temperature  result  of  25.2  °C  with  a  direct  contact  for 
the  same  heat  source,  it  can  be  deduced  that,  better  results  can  be  acquired  from  indirect 
contact  with  intervening  materials  with  higher  conductivity. 

The  faces  of  the  heat  dissipater  and  the  radiator  inside  the  CubeSat  were 
determined  as  bare  copper  because  of  the  low  IR  optic  properties  of  the  copper.  The 
reason  was  to  minimize  the  heat  flow  into  the  CubeSat.  On  the  other  hand,  some  of  the 
experiments  were  conducted  to  determine  the  effectiveness  of  surface  finish  of  the 
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outside  face  of  the  radiator.  As  expected,  the  results  showed  that,  optical  properties  of  the 
outside  surface  of  the  radiator  are  very  effective  for  the  design  success.  For  example, 
even  though  bare  copper  1U  size  radiator  resulted  in  56.1  °C  circuit  board  temperature  for 
a  test  configuration,  black  tape  covered  same  size  radiator  reduced  the  temperature 
46.5  °C.  Similarly,  covering  the  3U  size  radiator  with  black  tape  to  simulate  a  black 
painted  surface  reduced  temperature  to  26.4  °C,  from  the  bare  copper  radiator  result  of 
43.9  °C.  Even  though  the  black  photographic  tape  was  used  for  this  research  to  simulate 
black  paint  and  improve  the  optical  properties  of  the  radiator,  it  is  expected  that,  a 
polyurethane  black  paint  coated  surface  minimizes  the  contact  resistance  and  can  yield 
better  results. 

Radiator  size  is  another  important  factor  for  the  design.  ET-H  and  ET-J  results 
showed  that,  if  the  CubeSat  mission  allows,  the  radiator  size  could  be  increased,  such  as 
2U  or  3U,  to  acquire  better  results  for  higher  power  processors.  Even  though  it  was  not 
tested  in  this  research  because  of  the  time  constraints,  more  than  one  face  of  the  CubeSat 
can  be  used  for  radiators  and  the  heat  dissipater  can  be  connected  to  those  radiators  using 
thennal  straps  to  deal  with  very  high  heat  levels.  It  is  also  proved  in  the  research  that, 
using  a  deployable  radiator  is  another  alternative  in  case  of  high  heat  levels.  However,  it 
should  be  noted  that,  assuming  there  is  very  low  conductivity  through  deployable  radiator 
hinges,  increasing  the  conductivity  between  the  deployable  radiator  and  the  heat 
dissipater  box  is  essential  for  better  results.  If  relatively  long  straps  are  to  be  used 
between  the  heat  dissipater  box  and  the  deployable  radiator,  materials  with  very  high 
thennal  conductivity  may  be  selected  to  increase  the  perfonnance. 
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Finally,  orbit  simulation  predictions  proved  that  the  heat  dissipater  design  is  very 
successful  to  remove  the  heat  in  both  sample  cold  case  and  hot  case  scenarios.  The  design 
is  flexible  enough  to  let  the  designer  to  pick  a  radiator  size  and  number  of  straps,  which 
gives  the  best  results  for  the  situation.  In  case  of  very  high  heat  levels,  deployable 
radiators  are  also  preferable.  The  design  is  also  very  effective  when  a  direct  contact  is  not 
possible  between  the  high  power  processor  and  the  heat  dissipater  box.  Consequently, 
this  design  is  a  modular  approach  to  dissipate  heat  from  high  power  heat  sources  on  a 
circuit  board  in  a  CubeSat.  It  is  capable  to  prevent  excessive  temperature  condition  for 
the  heat-generating  component  and  to  minimize  heat  flow  into  the  satellite.  Due  to  its 
modularity,  this  design  can  be  used  in  any  1U  stack  in  a  CubeSat. 

Recommendations  for  Future  Research 

The  design  could  not  be  tested  with  various  materials.  Using  materials  with  higher 
thennal  conductivity,  such  as  OFHC  copper,  instead  of  commercial  copper  plates  in  the 
box  and  radiator  design  would  be  a  better  option.  Additionally,  straps  with  higher  thennal 
performance  than  copper  ground  wire  braids  would  improve  the  affectivity  of  the  design. 

All  experiments  were  conducted  with  1  mm  thick  copper  radiators  in  this 
research.  Thicker  radiator  perfonnance,  especially  for  larger  scale  radiators,  2U  and  3U, 
and  surface  finish  effect  of  the  radiator  are  also  important  factors  should  be  examined. 

Thermal  strap  connections,  such  as  bolting  the  straps  using  a  highly  conductive 
intervening  material  or  using  thermal  paste,  are  another  aspect  would  be  searched.  Using 
radiators  on  more  than  one  CubeSat  face  to  deal  with  very  high  temperatures  would 
provide  greater  confidence  to  design.  However,  the  thennal  performance  of  relatively 
long  straps  would  be  another  challenge  in  this  configuration. 
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The  heat  dissipater  design  performance  was  determined  only  by  using  two 
different  circuit  boards  in  a  sample  6U  CubeSat.  Experiments  with  various  payloads  in 
different  size  CubeSats  would  assure  the  modularity  of  the  design. 

Finally,  since  a  few  computational  simulations  were  conducted  only  to  determine 
the  mechanical  strength  of  the  heat  dissipater  design,  it  should  undergo  detailed 
computational  mechanical  analysis  and  physical  testing  to  determine  the  flight  acceptance 
level  and  the  ability  to  operate  on  orbit. 
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Appendix  A:  Technical  Drawing  of  the  Heat  Dissipater  Box 

Technical  drawing  of  the  heat  dissipater  box  is  shown  in  Figure  47.  All 
dimensions  are  in  inches.  Material  of  the  faces  (1  and  4)  are  copper  and  the  edges  (2  and 
3)  are  aluminum. 
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(a)  Exploded  View  of  the  Heat  Dissipater  Box 
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(b)  The  Top  Face  of  the  Heat  Dissipater  Box 
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(c)  The  Bottom  Face  of  the  Fleat  Dissipater  Box 
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(e)  The  Second  Edge  of  the  Heat  Dissipater  Box 

Figure  47:  Technical  Drawing  of  the  Heat  Dissipater  Box 
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